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FOREWORD 


The  computations  described  in  this  report  were  performed  in  the  Fire  Control  Formulation 
Branch  (K41),  Submarine-Launched  Ballistic  Missile  (SLBM)  Research  and  Analysis  Division, 
Strategic  and  Strike  Systems  Department.  This  report  is  a  compilation  of  all  the  important 
mathematical  results  and  their  derivations  pertaining  to  the  application  of  the  Lambert  Problem 
to  inverse-square  gravity.  In  particular,  expressions  for  the  correlated  velocity  and  the  null-miss 
vector  in  terms  of  known  parameters  are  derived.  The  treatment  is  carried  out  in  detail  for 
elliptical,  parabolic,  and  hyperbolic  trajectories.  The  appendices  include  background  material 
needed  for  the  main  report,  alternative  derivations  of  or  expressions  for  some  key  variables,  and 
a  derivation  of  an  expression  for  the  rate  of  change  of  the  null-miss  vector. 

It  is  intended  that  this  report  will  serve  as  a  comprehensive  reference  guide  for  all  those 
who  desire  easy  access  to  the  Keplerian  trajectory  equations  simulating  the  trajectory  of  a  free- 
fall  reentry  vehicle  in  an  inverse-square  gravitational  field. 
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INTRODUCTION 


The  correlated  velocity,  vc ,  is  that  velocity  which  will  carry  a  reentry  body  on  a  free-fall 
trajectory  from  the  initial  position  (at  release),  ij ,  to  the  target  position,  r2 ,  in  a  specified  time, 
tf .  A  guidance  computer  steers  the  missile  in  a  direction  which  will  cause  the  difference 

between  the  correlated  velocity  and  the  true  velocity  to  approach  zero.  When  this  difference, 
which  is  called  the  velocity  to  be  gained,  is  zero,  the  reentry  body  is  released. 

The  Lambert  Problem  is  to  determine  vc ,  given  r, ,  r2  ,  and  t  j .  It  is  well-known  that  there 
is  no  closed-form  expression  for  vc  in  terms  of  the  given  parameters.  However,  vc  can  be 
determined  by  employing  an  iterative  procedure  that  utilizes  the  Newton-Raphson  method. 

In  this  report,  three  candidates  for  the  iteration  variable  will  be  considered.  They  are: 

AE 

(i)  v2,  where  x  =  ,  AE  being  the  change  in  the  eccentric  anomaly; 

(ii)  vcd ,  which  is  the  tangential  component  of  vc  ; 

h 2  .  .  . 

(iii)  p,  known  as  the  semi-latus  rectum,  where  p  =  — ,  h=  r,xv  ,  and  p  =  GM,  G 

M 

being  the  universal  gravitational  constant  and  M  the  mass  of  the  earth. 

It  will  be  shown  that  x,  vc9 ,  and  p  are  related  by  the  equations 


2  2 
ri  vce 


2  6 R 


2 r,  p  sin  — 
2 


r\  +  r2  ~ 


9 


qr2  cos  x  cos 


where  0R  is  the  range  angle  defined  by 


0R  =  arccos^— =  arccosl^  •  r2),  0  <  0R  <  n . 
r\r2 

Using  the  above  relations,  we  will  derive  the  equations  for  an  elliptical  trajectory  in  terms 
of  each  of  the  above  iteration  variables.  We  will  then  deduce  the  equations  for  hyperbolic  and 
parabolic  trajectories,  and  demonstrate  how  the  type  of  trajectory  can  be  determined  from  the 
given  parameters.  Since  we  will  be  using  results  that  are  derived  in  Appendix  A,  it  is  advisable 
to  read  Appendix  A  before  delving  into  the  derivations  that  follow.  Note  that  all  the  equations  in 
this  report  are  based  on  the  hypothesis  of  an  inverse- square  gravitational  field. 
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TIME  OF  FLIGHT,  CORRELATED  VELOCITY, 
VELOCITY  AT  IMPACT 


It  is  shown  in  Appendix  A  that  the  time  of  flight  for  an  elliptical  trajectory  is  given  by 


tfe  =  Y {E2  ~  E'  ~  e  (Sil1  E2  ~  ^  )}  = 


r 


A E  E,+E- 


A E  -  2e  sin - cos  — 1 - - 

v  2  2  , 


=  2j—  {x-esinxcos 


(*+£i)}>  x=^~- 


(i) 


Now 


+  r2=a(  l-e  cos  El )  +  a  (l  -  e  cos  E2)  =  a{  2  —  e  (cos  £j  +  cos  E2 )}  =  2a  {l  -  e  cos  x  cos  (*  +  £>)} 

(2) 


;cos(x  +  £j)  =  |  1  -  r‘  +  Vl 
' '  ''  2a  j 


secx . 


,R  9n  —0,  6,  0 t,  .  6,  .  02 

cos—  =  cos  — - - L  =  cos  — -  cos  — —  +  sin  — -  sin  — m 


2 


2  2  2  2 


a  ( 

l-e)  £j 

COS 

_  la(l-e)  E,  MH  £2  a(l  +  e)  £j  a(l  +  e)  E 


2  V  n 


■  cos  — —  + 


sin- 


2  V  r, 


-sin- 


ci  \  / 1  \  .E,  Ej  / 1  \  .  E,  .  Ej 

(l-e)  cos — cos  — -  +  1  +  e )  sin  —  sin  — - 
r,K  l  2  2  y  2  2 


12 

fl  C  E2~E1  E\+El\  a  (  (  ,  77  11 

cos  — - —  e  cos  — - -  >  =  , —  |cosx-gcos(x  +  Ejj. 

2  2  J  Jrtr2 


■■■  n 


Q  Y 

+  f2  -2-^/ip^cosxcos^-  =  2a  sin2  x  =>  a  =  — csc2x 


0 , 


/•j  +  r2  -  2^/r,r2  cos x cos  — 

2  y 


Q 

For  the  sake  of  brevity,  put  S  =  rl+r2  -2 ^rxr2  cos  x cos  as  this  expression  will  be  used 

frequently.  Then 


a  =  —  S  esc2  x . 
2 


(3) 
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The  iteration  continues  until  tje  is  infinitesimally  close  to  l  j ,  at  which  point  we  obtain  the 
correct  value  of  x .  Usually,  four  iterations  are  sufficient  for  convergence. 
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Proposition:  The  correlated  velocity  is  given  by 


v  =hzA 

C 

8 


where 


T  o  T  "R 

f= — ^  +  2  —  cos  v  cos — — ,  g 


0,, 


ri  V  ri 


2  r.r^S  0R 

— u^~  COS  — 

JU  2 


Proof:  Since  rj ,  r2 ,  and  vc  are  coplanar  vectors,  each  of  them  may  be  written  as  a  linear 
combination  of  the  other  two.  Let  f2  =  /r,  +  g\c .  Then 


Let 


r2  x  vc  =  r2  xr,  =  /  r,  x  vf  =  fh,  |r,  xr2  =  g  xvc  =  gh . 


*i 

r  = 

Ti 

>  ri  = 

yi 

0 

0 

with  similar  expressions  for  r2  and  r2 .  Then 


Jtj  =  rx  cos  0X  =  a  (cos  Ex  -  e) ,  y,  =  rx  sin  6x-b  sin  Ex , 


h 

with  similar  expressions  for  x0  and  y2  .  Using  the  fact  that  E  =  — ,  we  obtain 

br 

xx  =  -ciEx  sin  Ex=~— sin  Ex ,  yx  =  bEx  cos  Ex  =  —  cos  Ex , 
brx  rx 

with  similar  expressions  for  x2  and  y2  .  Hence, 


/  =  — — -  =  — ^2-  =  —[cos  Ex  (cos  E^-e)  +  sin  Ex  sin  E 2 } 

h  hr. 


=  —[cos  (E2  ~Ex)  -ecos£’,}  =  —  <  cos  (E2  -Ex)-\  +  — 


1-— [l-cos(£’2  -Ex )}  =  !-— sin2  jc  =  1-  —  =  -  —  +  2  ^  cos  x  cos  — , 
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g  =  — - —  =  =  —{(cos  £j  -  e )  sin  E2  -  (cos  E2  -  e )  sin  £j } 

h  hh 


3  I  3 

Id  r  „  d 


—  {sin  (E2  -  E 1 )  -  e  (sin  E2  -  sin  £j )}  =  2  —  sin  .*  {cos  v  -  ecos  (x  +  E1 )} 

M  \  M 


f\r 2  0K  ,,  ar.n  .  0  2 r.r2S  0,. 

“ —  —  -  0  — —  sin.rcos  —  =  — — —  cos  — — 


=  2  J — sinx  — — —cos  — —  =  2 

jU  a  2  V  jU 


M 


The  magnitude  of  the  correlated  velocity  is  given  by 


f  2  l) 

r 

r2  o 

vc=,M 

— 

=  J 

- — r 

V 

l«U 

where  a  is  given  by  Equation  (3). 

Proposition:  The  velocity  at  impact  is  given  by 

v,  =  /r  +  gx 

where 

cos  x 


/  __zr^Lj2fiS,  g=-  —  +  24  —  cos x cos . 


fV~2 


r2  \  r2 


First  proof: 

|v,  x  vc |  _  \*2  x  ^  |  _  x2y1  -  xxy2  _  ah 


f 


h  h  h 

•  |fi  xvj  kxfj 


sin(£’2  -£j)  =  -  sin  (2,v)  =  _C0S  A 


rir2 


rir2 


2a  .  2  r,  „  I  r,  0D 


h 


h 


1-— sinz  x  =  1-  — = — -  +  2  j— cos.vcos  — . 

2 


r2 


r2  r2  V  r2 


Second  proof:  By  considering  the  reverse  trajectory,  we  have 

r,  -  /?2 


v,  =  — 


where 


,,  ,  S  2r^S  0R 

f  =1 - ,  8  =  -^^cos  ~^~  =  g  ■ 

V  E  2 
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Now 


v,  =  M  +  gv,  =  M+-(r2-fi,Y 

8 


f- 


fg_ 

8 


•  \ 


ii+if2  1  s-f 


8 


r  ’ 


8  8  8  8 


The  second  equation  yields  immediately  g  =  f.  The  first  equation  yields 


/ 


fg- 

8  8  8 

C0SXj2 JS. 


f 

f  <0 

r  s) 

1 

1-  — 

1-— 

-1 

l 

l  rJ 

V  r2> 

J 

8^2 


■{S-rl~r2)  =  — 


IS 


cos v cos - 


0, 


<?V?ir2 


Third  proof:  The  expressions  for  /  and  g  may  also  be  obtained  from 

3/  a  •  ^8  h 


f  dE2  El'  8  dE2  ~2’ 


where 


2a  .  2  E2-  E  a 3  r  . 

/  =  1 - sin"  — - 1  1 — 


8  =  J~  (sin  (E2  ~  Ei )  -  e  (sin  E2  -  sin  Ex )},  E2  =  —  J—  . 


2  V  E 

The  magnitude  of  the  velocity  at  impact  is  given  by 


vt=jE 


r  2  1A 


Vr2  aJ 


r  2__J_A 

v  r2  h ) 


NULL-MISS  DIRECTION 


The  null-miss  direction  is  the  direction  along  which  a  velocity  error,  Avc ,  at  release  can 

exist  and  the  reentry  body  will  still  hit  a  target  on  a  rotating  earth.  It  is  assumed  that  the 
atmosphere  is  moving  with  the  earth  and  no  winds  exist.  The  null-miss  direction  is  defined  by 
the  unit  vector  PI)K  ,  where 


Av^ 

At, 


ldk 


" 

J 

A  vc 

_  c 

“  |Avc| 

Mf 

=  -l. 
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Since  tfe  =  l  fe  (x.  0R ) ,  wc  have 


dtfe  dtfe 
Atf  =—Ax  +  z^rA0R, 


dx  dO, 


R 


where 


ri  xr2 


A 0R=£l-  nA tf  =  -£l  n,  n  =  ^ ^  =  r,  x  r2  esc  0R , 


ri  xr2. 


and  fl  is  the  vector  representing  the  earth  rotation  around  itself. 


At, 


dt 


fe 


Ax  =  -  d°R 


A  0R  1  +  ^A  0 
R  dG B  ‘ 


dt 


fe 


dt 


fe 


dx 


dx 


dt 


where  — —  is  given  by  Equation  (5),  and 
dx 


dtf  1  /'r,  .  0 

_ l _  — _ _ L_± _  c  i  ‘ 


sin- 


dOR  2  V  2 JUS  2 


(ri  +  i'2){3csc2v(l-.x;cot.x;)-l} 

Q 

“6 ^Jr\r2  cos  esc  x  cot  x  ( 1  —  x  cot  x) 


_  r  —  fr 

Similarly,  by  taking  differentials  of  both  sides  of  the  equation  vc  =  — - - ,  we  obtain 

g 

Ar0  A/  _  Ag  _ 

Avc  =  -  —  Tj - ^vc, 

8  g  g 


where 


Now 


Ar,  =  Qxr7Atf  =  -ftxr„  A/  =  ^-Ax  +  ^—A0R,  Ag  =  ^-Ax  +  ^-A0K. 
2  2  7  2  dx  d0R  R  dx  d0R  R 


df  jr\  .  0R  df 
-f—  =  -2  —  sinvcos  — ,  A— 
dx  \j  /,  2  d0R 


y\  .  \s  n 

-  —  cos  jc  sin- 


?r 

2 


dg 

dx 


2  .  2  @R 

-r^sinvcos 


dg_  A 


AS 


2  ’  d0„ 


sin 
2  juS  2 


0 B 


rtr2  cos  v cos  -  r{  -  r2 


(6) 
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Hence, 


f 


ov 


AS  =  J— 


2 

juS 


A/  =  -  24P 

>i  v 

r,r2 sinvcos2  —  Ax  +  1  '1'2  c';" 


0„  A0P 


sin  x  cos  —  Ax  +  cos  x  sin  —  — - 
2  2  2  , 


r,n,  .  0D 
■sin 
2jlS  2 


0R 

r\r2  cos  x  cos  -  rt  -  r2 


A0„ 


2  •  2  Or  K 

- r.n  sin  *  cos'  —  Ax  + 

A  g_i>uS12  2 


r.r2  .  0R  r 
1  -  sin  — 


\j2jLiS 


V 


,  0R 

r\r2  cos  x cos  -  rx  -  r2 


a  oD 


g 


2  r,r,S  0R 

I  M  2 


Or  a  1  0/ 

- sin  x  cos  —  Ax  H - tan  — 

S  2  2S  2 


0R 

rYr2  cos x cos —  rt  —  r2 


A  Or, 


Or  -Or  A$.  ' 

sin  .x  cos  —  Ax  +  cos  x  sin — — - 

2  2  2  j 


Or  A  Or 

-tan-  R  R 


2  2 


.  .  Or  A  -On  A  On 

Putting  z  =  sin  .x  cos  — -  Ax  +  cos  .x  sin  — - —  ,  we  obtain 


2  2 


Av„ 


flxr,  2 z  pT- 
- -  +  —  — r,  + 

8  g  V  ri 


Or  A  Or 
tan-  R  R 


2  2 


V.  . 


(7) 


Equation  (7)  may  be  written  in  terms  of  r, ,  0, ,  and  n  which  form  a  right-handed  orthogonal 
coordinate  system.  We  have 


Now 


Avt  =  ( AV  f, )  ?,  d 9, )  6,  +  ( Avr  •  A)  n  . 


Qxr,  r,  =  Q  ■  r2  x r.  =  -r,Q  •  ?,  x r2  =  -r,Q  ■  n sin  0R  =  r,  sin  0RA0R , 


ilxr,  -0,  =fi  r2x01  =  fir2x^ — E cos 0R  _  r,  x r2 cot 0R  =r2fl-ncos0#  =  -r2cosORAOR. 

sin  0L, 


Qxr,  =  (iixr2  - f , ) r,  +  (0xr2  - 0, ) 0,  +  (fixr,  -n)n 

=  r2  sin  0RA0Rrx  -  r2  cos  ORA0RQl  +  (0xr2-n)ii. 
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-  _  r2  ~  fii  _  r2  (e cos  ^  +  9,  sin  ^ )  -  frxrx  _  ( h  COs  0 R  -  frx )  r ,  +  a  sin  0^9, 


v„  = 


g  g  g 

f  r -  6K  V 

r2  cos  0R+r2-  2yjt\r2  cos  x  cos  —  r,  +  r2  sin  0fl0, 

2  y) 


'  0 


_  v 


A  2  ^  ff  A  /  A  «  ■  ^p  A 

2a  cos  —  ~2^r\r2  cos x cos -y-  Jr,  +  2r2  sin-^-cos-^Bj 


2k  aS1  0r 
1  -  cos  — 


A  0fi 
—  cos— -cos  X 

vJri  2 


T+. 


2/^r2 


sin  — 0,. 
^  1 


rtS  2 


Substituting  the  expressions  for  0xr2  and  vc  into  Equation  (7),  we  obtain 

1  2^ 

Avc  = - jr2  sin  0KA0Kr]  -  r2  cos  8RA0RQX  +  (O  x  r2  •  n)  ii}  +  -  ^  1  ~ 


g 


+ 


0R  AO..  Jrr 
tan-2- — 


2  2 


vVri 


h, 

2 


K.  \y  n 

—  COS  —  -cos  X 


,  2jur,  .  0R  ~ 
r,  +  ~  sin  — 0. 

V  t]S  2  1 


,2  9r 


Ax 

-  ^  ,  e ^ 

rx  sin  x  +  r2-  r2  cos  x  cos 

v  "  2  J 


A0„ 


(8) 


+ 


li 


in  OrAx+  T%<|(?i  +  r2)cos%-7^ 


2S 


A  sin  x  sin 


3 '2 


f 


'  2rtS 


cosx 


0, 


\ 


1  +  cos2^ 

V  2; 


A0,„ 


—  (Ox  r,  •  ii)n. 

gV  ' 


REENTRY  FLIGHT  PATH  ANGLE,  ANGLE  OF 
ATTACK,  AND  YA 

These  three  angles  are  computed  at  reentry,  i.e.,  at  400,000  ft.  Let  r2  be  the  position  vector 
corresponding  to  that  altitude.  Then 
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l- H 

E'2  =  2^-arccos - — ,  x'  =  —  (E2  -£j),  f'  =—  {2.r/-<?(sin£,j  -sin£j)}, 

el  h 

/'  =  1  sin2x',  ,g/  =  ^{sin(2.x:/)-e(sin£2 -sin£,)},  r'  =  /f,  +  g' vc, 


>i 


/'■  /  (I h  .  f  ~  1  2fl  .  2  /  —  /  .  ■/—  -V  /a.  ~/\ 

/  =  — : — ;sm(2-*  ),  8  =1-  —  sin  X,  V,  =  /r,  +  gvc,  =  arccos(r, -r2). 


The  velocity  at  reentry  takes  into  account  the  earth  rotation.  It  is  given  by 

v;  =  v;-nxr2\ 

The  reentry  flight  path  angle  is  the  angle  between  x'r  and  the  plane  normal  to  the 
instantaneous  local  geodetic  vertical,  h  (Reference  2).  It  is  given  by 

71  (C 

Yr  =  — —  arccos(h  •  v  J. 

The  null-miss  angle  of  attack  is  the  angle  between  PDK  at  release  and  v' .  It  is  given  by 

a  =  arccos(PDA.  •  v' ) . 

If  vr  is  the  analog  of  v'  at  impact,  then 

vr  =  v,  -flxr. 


and  the  zero-angle  of  attack  is  given  by 

a0  =  arccos(vr  •  v' ). 


Let  (n,,u,vy)  be  a  right-handed  orthogonal  reference  frame,  where 


/V  / 

n 


k'xv'  . 
—  /  —  /  ’ 
r2  X  V  r 


then  u  =  v'xn.  Let  YA  be  the  angle  between  u  and  the  projection  of  PDK  on  the  nu  plane,  to 
be  measured  clockwise  from  u  in  the  n  u  plane.  It  can  easily  be  shown  that 


YA  =  arctan 


P DK  ' 
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TIME  OF  FLIGHT  IN  TERMS  OF  v 


C0 


From  the  expression  for  vc  on  page  gj  we  have 


vce 


P^sin^ 
I  rtS  2 


.  i  9r 

Q  2 ur2  sin“  — 

c  ,  0  / -  Or  9 

S  =  rx  +  r2  -  2 •x/r1r2  cos  x  cos  — -  = 


r,v 


1  V  C0 


(9) 


ri  ( ri  +  r2 )  -  2FF  sin 


•  2  0R 


cos  v  = 


2  0R 

'\\lr\r2vce  cos  2 


(10) 


sinx  =  vt  —  cos”  x 


j4/i  Vi  cos2  -  j#i  (#;  +  r, )  v2e  -  2 fir2  sin2  ^ 

-  I -  2  0R 

2riVrir2ve0cos  2 


2  4 


(p  +  -  2rtr2  cos 0* )  +  4//rtr2  (#;  +  r2)  v2*  sin2  ^  - 4y£Tr22  sin4  ^ 

2riv/Vicosy 


_v 


2/^2  sin2  y-/pp 


'i  +  '2  - 


cos- 


0 , 


0D 


Q 

n  +  h  +  2 Jr,r2  cos -p  |  -  2/ir2  sin2  -f- 1 


2  0 R 


'\yr\r2Vce  cos 


Put 


F  =  l2//r2sin  y-rp; 


■  2  0R  ^  2  C 
" '  \:0 


-  / - 

ri  +  r2  -  2v rir2  cos  y 


/-  +  r2  +  2pp  cos  ^  -  2/ir2  sin2  ^  j . 


Then 


sin  x  = 


Vf 


2  ^  ’ 
'\^Jr\r2Vt;0  COSy 


(11) 
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Observing  the  expression  for  F,  we  see  immediately  that  the  above  inequality  implies  that  F  is 
positive  for  elliptical  trajectories. 
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vc6  is  updated  by  the  Newton-Raphson  method,  where  a  good  initial  estimate,  according  to 
Reference  3,  is 


vce  = 


O  -2  Or 

2/ir2sin 

r2  sin  0R 

+ 

rl 

(  Q  \ 

n  +  r2+2ylr^cos-£- 

tf 

which,  by  the  above  inequality,  holds  provided  that 


1  q2  +  f~2  ~~  2qr2  cos  0R  ^ 

I  f  ^ 


lf 


2fi 


rl+r2-  2^jr{r2  cos +  J t \  +  r2+  2 ^r\r2  cos  ^ 


vcr ,  the  radial  component  of  vc ,  may  be  expressed  in  terms  of  vc9 .  From  the  expression 
for  vc  on  page  lo]  we  have 


v„  = 


lr2  0R 

\—  cos  —  -cos  x 

\  V  ri  2  J 


Substituting  Equations  (9)  and  (10)  into  the  above  expression,  we  obtain 

r  v 


cot  0R  — -  esc  0R 


v 


j 


u  0R 

vc9  ~l  “  tan  . 
r\Vc$  2 


Hence,  vc  =  v^r,  +  ve001 ,  where  =  nxq  =  r2  cscdR  -fj  cotdR  .  By  considering  the 
reverse  trajectory  and  using  the  fact  that  h  =  rlvc9  =  r2vt9 ,  we  obtain  for  the  velocity  at  impact 


where 


v  =v,E+v,„0,, 


ri 

v,e=-vce. 


f 

> 

0R 

tan  —  = 

f 

—esc  0R 

1 

0 

0 

r-f 

** 

CSC0R 

1 

I-" 

0 

0 

r-f- 

vte 

u 

9 

W* 

2 

K 

r2  ) 

u  0R 

-2— tan  — , 

r,V„n  2 


02  =  n  x  ?2  =  n  x  (r,  cos  0R  +  0,  sin  0R  j  =  0,  cos  0R  -  r,  sin  0R  =  r2  cot  0R  -  r,  esc  0R . 
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We  have 


TIME  OF  FLIGHT  IN  TERMS  OF  p 


h  _  Jjup 


ri  r, 


Substituting  this  value  of  vc6  into  Equations  (13)  and  (14),  we  obtain 


lfe 


\P  rir2  Sin^jg 
F' 


2 rxr2  {rl  +  r2)  sin2  y-/?  (r;  +  r22  -  2rxr2  cos 0R ) 


2 pr,  n  sin"  6,, 

+  1  — - -  arctan  - 

VF 


(rt  +  r2)  p  —  2  rtr2  sin 


2 


(15) 


dtfe  dtfe  dvc0  1  U l  dt 


fe 


dp  dvc0  dp  2/,  ]]  p  dv 


c8 


=  y, |  P2  (K  +  ri~  V2  cos 0R )  -  4ifr22  sin4  ^ 


*  fe  p  2rp  sin3  0,, 

1 -+.  - - 1  2  , —  R  arctan  - 


Vf7 


2p  Vju  fVf7 


(/,  +r2)p  —  2 1] r2  sin 


•  2  6r 

2  J 


8/|4r,4  sin3  0R  sin2  — 
IP  12  *  2 

77 '2 


(16) 


where 


77'  j-)  •  2  6r  f 

F  =  j  2rtr2  sin  y  -  p 


rx  +  r2  -  ZJr \r2  cos- 


0 B 


v 


+  r2  +  2  Jr,r2  cos 


—  -2rn  sin 

2  1  z 


■  2 


—  F 
F2 


2  2 
V  V 

p  is  updated  by  the  Newton-Raphson  method,  and  a  good  initial  estimate  is  p  =  — — — 

M 


where  vc9  is  initialized  as  on  page  [j_3 
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EQUATIONS  FOR  A  HYPERBOLIC  TRAJECTORY 


From  Appendix  A,  we  have 


„  <?  +  cos  0  .  ^  yjl-e2  sin# 

cosE  = - ,  sinU 


1  +  ecos# 


1  +  ecosO 


From  the  expression  for  sin  E.  we  observe  that  since  e  >  1  for  a  hyperbola,  E  is  imaginary  in 

AE 

this  case,  i.e.,  E2  is  negative.  Since  x  =  -y- ,  it  follows  that  x  is  imaginary,  and  hence  x2  is 


negative.  Also,  a 


h2 


<0. 


I u(l~e 2) 

Let  H  =  \!~E2 ,  y  =  yj—x 2  ;  then  E  =  iH,x  =  iy,  sin  E  =  /sinh  H ,  cos  x  =  coshy, 

sin  x  =  i sinh  y,  tan  x  =  i  tanh  y,  xcsc x  =  y  csch  y,  xcot  x  =  y  coth  y,  esc2  x  =  -csch2y  ;  also, 

arctan(zw)  arctanh  u  1  ,  1  +  u  „  i  i  ,  TT  .  ,  „ 

v  -  for  \u\  <  1.  Using  the  fact  that 


IU 


-In- 

2  u  1  —  u 
2 fir2  sin 


2  @R 


2jiir2  sin 


2  &R 


VC0 


> 


rx  +  r2  -  2yjrlr2  cosh  y  cos  ^ 


e A  / 

r, 


n  +e 


r{r2  cos 


we  see  immediately  that  F  and  F'  are  negative.  Putting 

4—f 


u 


Q 

i  ( ■ ri  +  F )  vle  -  2AL  sin2  y  ( ri  +  ri )  P  ~  2r7i  sin  2 


2 


Q 

and  noting  that  rx  (rt  +  r2 )  v20  -  2 // r2  sin2  -^->0,  it  is  easy  to  verify  that  u  <  1.  We  have 


-F  ~  j/i  U  +  r2)v2e  - 2 iir2  sin2  ^  j  =  -4r,V2v^  cos2  ^  <  0 

r 


J  +  r,  ( rt  +  r2 )  v20  -  2 /ir2  sin2  j  1 F  - 1]  ( rx  +  r2 )  v2cd  +  2 fjr2  sin2  y  j  <  0 

=,V^<l-|(l-1  +  r2)v2  - 2^2sin2|  =*«<!. 
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_  Q 

Putting  S'  =  r{  +  r2  -  2-s]rlr2  cosh  ycos-^-,  we  obtain  the  following  expressions  for  the 
hyperbolic  time  of  flight: 


tfli  =a^—{H2  -  Ht  -<?(sinh//2  -sinh//,)} 

ry  r  _  q  j 

=  I— csch2  y  |(r,  +  r2 ) (cosh  y  —  y  csch  y)  +  2 cos -^-(y  coth  y  - 1) j 


rir2vC0  sin  0R 


_  /  \  .  2  ^  Tt  9  /  2  2  r*  ri  \  ^ sin"  On  1  H ~  li 

2/^2  ( rx  +  r2 )  sin"  -f  -  pv;,  r;  +r2-  2  rxr2  cos  0R  +  -  f= - In  - - 

2  J-F  1 -u 


4^ 


Ip  r,r2sinOR  „  ,  \  ■  ^  0R  ,  ,  2  „  _  \  pr  V,,2  sin2  ,  1  +  w 

= - — — - — M  2rr,  (r,  +  r2 ) sin"  —  -  p[  r:  +  r"  - 2r,r~  cos +  1  ~ — = — -In - 

]IjU  F'  1  12Vl  27  2  V  1  2  12  1  —  m 


Also, 


— —  =  ,  i  csch  y 

42JitS' 


(/,  +  r2)2{3csch2y(ycothy-l)-l} 


+2r,r2  cos^  —  {3y  coth  y  (2  csch2  y  + 1)  -  6  csch' y  -  5j 


+-J44  (ri  +  ri ) cos 12 csch2 y  +  l)cosh  y  -  3y  csch  y  (4 csch2 y  +  3)j 


=  j  (r,2  +  r22  -  2 rxr2  cos  )  -  4/r  r22  sin4  ^ 


tfli  _|_  2/// 1  r2  vr6)  sin  0R  ^  1  +  « 


V  vce 


f4^f 


l  —  u 


On 


dtfh  1 


=  \  P2  (4  +  r2  -  2rir2  cos  0R )  -  4r,2r22  sin4 


dp  F 


lfh 


16//  r[r2 vcd  sin3  0R  sin' 

T2 

p  pr2  sin3  0R  1  +  u 


+ 


2  p  Va  f'4^f'  i -«y 

8r,4r4  sin3  sin2  — 
IP  12  2 

77/2  ' 


at 


a#,. 


_ 


1  / ,  r,  .  0,, 
— ,  ",  sin  — 
2^2/pS'  2 


(ri  +  p){3csch2y(ycothy-l)-l} 

Q 

-69 J~t\4  cos  csch  y  coth  y  ( y  coth  y  - 1) 
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y2,  vfe,and  p  may  be  updated  by  the  Newton-Raphson  method  as  follows: 


2  _  2,.2y. {>,-<„ (?■)}  ,  > 

y»+ 1  ^  \  >n+\  v  Afl  /„  ^ 

fh 


dt 


dy 


f  dt  ^ 

Ulfh 


n+ 1  r  n 


Pn+- 


rdt  ^ 

Ulfh 


\  d VC0  Jvce=(vce) 


dp 


where  the  usual  initial  guesses  may  be  used. 

The  correlated  velocity  and  velocity  at  impact  are  given  by 


vc=^,  v,  =  fhr1  +  ghyc , 

Sh 


where 


/*=-— +  2.  —  coshycos-^,  gh  =  fh  =  gh  = 

>]  V  >\  2  V  A  2  rir2 


r.  _  r. 

-—  +  2  —  coshycos  — . 


Alternatively,  vc  =  +  vce0P  v,  =  vtrr2  +  vr0Q2 ,  where  vcr,  vc0,  vlr,  and  vte  are  given  on 

page  13.  Also, 


vc  =  AP 


f  2  1A 


v  r\  a  j 


f  2  1  A 

— +  n 

Vri  \a\j 


’  Vr=JP 


^  2  1A 


V  r2  aJ 


1 


A 


f  2  1  A 

— +n 

\r2  M  J 


The  null-miss  vector  is  given  by 


Av  =_^xF2^2z/  Vit 


Si, 


-  H - .  I— r,  + 

Sh  v  'i 


tan^A^_V^ 

2  2  5' 


2p  .  (  .  ->  I -  .  ^1  A 

— ^  sinh  y  rx  +  r2  sin  —  ~yjrxr2  cosh  y cos  —  Ay 
S  V  2  2  7 


up  .  0,, 

+  - _  sin  — 

\  2r1S  3  2 


-r{  sinh2  y  +  r2-  ^r\r2  cosh  y  cos 


eD 


a  e„ 


+ 


yy'i  sinh  y  sin  6RAy  +  j  (r,  +  r2 )  cos  ^  cosh  y 


f.  2  0R) 

l  2  J 


A  0,, 


where 


—  (ilxr2  ■  n)n, 

Si, 


A 0R  =  -Q.n,  Ay  = 


1  +  ^ALA6R 
d Or 

dtfh 

dv 


/  .  0R  .  Or  A  6r 

z  = -sinh y cos  — Av  +  cosh y sin— — — 
2  "22 
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EQUATIONS  FOR  A  PARABOLIC  TRAJECTORY 


Since  the  time  of  flight  is  a  well-behaved  function  of  x2,  it  follows  from  Equation  (4)  that 
the  time  of  flight  for  a  parabolic  trajectory  is  given  by 

tfp  =limt/f  =lim  esc2  *  jOl  +  r2)(xcscx-cosjt)  +  2.N/^cos^L(l-Jccotx)j 

fl  i  2  V  f  ( 4V  /3/^24  Vl 

hm  -rt - 1 - h ...  r,+r,)  —x - I-...  +2Jr.r,  cos— - 1 - I-...  > 

lx2  3  15  JV1  2 \3  45  J  Vl2  2  l  3  45  J 


1  l2S0  f  I —  dR  3  / —  0R 

=  -J—  r1  +  r2  +  y/r1r2cos-£-  ,  where  S0  =  rl  +  r2-2y]rlr2cos-f- . 

•'ll  V  2  )  L 


Also,  from  Equations  (5)  and  (6), 


dt*  1  f,  \  2  1 0R  I — /  ,  0R  ] 

lim— JT  =  — ,  <( r.+r,)  -r.r7  cos  —  +  Jr.r2  ( r  +  r,  )cos—  >, 

~°3*2  5^2/^  l  2  V12V.  2t  2  J’ 


..  or  32  2S0 

lim— —  =  lim#-  = - 2 — 

x~>°  dvc8  -M°  dvcg  5 jUr2  sin  dR 

dx2 


in  +  h)2  -  Vi  c°s2  y  +  (#;  +  r2)cos^ 


y1  =  iim- = - - - a-J^—  I  (ri  +  -  rir2 cos2  y  +  v/^E  (>i  +  e)cos^4, 

^  ~°  ^  5rtr2  sin  sin  ^  V  l  2  2  J 


-^deR  2^2^ 


The  correlated  velocity  and  velocity  at  impact  are  given  by 


Vf  =  *2  ^  .  v,  =  /0*i  +  g„vc , 

£o 


where 


,  S0  12  r^S0  eR  ■  j2//S0  .  s0 

=  1-  — .  g0  =  J-LLJLcos1*-,f0=-l-^,  g0  =  l— 2-. 


h 

Since  r/  =  — 7 - —  and  e  =  1  for  a  parabola,  it  follows  that  a  =  °°  for  a  parabola.  Hence, 

Ml-*  ) 

vc  =  lim  ju  =  —  ,  which  is  known  as  the  escape  velocity.  Similarly,  vt  =  —  . 

a^°°]  Pi  aj  V  ri  V  r2 
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The  null-miss  vector  is  given  by 


A_  fixr,  2 z0  jr2  _ 
Av„  = - -  +  — 2- .  /— r,  + 


So  So 


tan- 


0K  A0R 
2  2 


J 


where 


'2  S0\ 


/,  +  r2  sin 


•  2  0, 


2  VR 


9r 

r.r0  cos  — 

I  Z 


A.r2 


///;  .  9r 

-  sin  — 

'2/A3  2 


/,r2  cos- 


0* 


v 


Atf„ 


+ 


—  -^-r,  sin  0D  Ar2  +  '  ^ 
2\2Sl 


t  \  dR  I - f 

,  Ar.  +  r2 ) cos  —  —  Jr,r~ 

2^1  2  Vl2 


6, 


1  +  cos2  — 

v  2  J 


AO, 


—  (Hxr,-n)n, 

So  V 


Ar"  =  -- 


dtf. 

1  +  lim — —A9r 
?)0K  R 

r  •  ^  fe 

lim  ' 

-v^°  dx 


1  f 

Z’’  =  2V 


’  ^0 


cos  — Ar"  +sin  —  A6R 
2  2 


The  equation  for  the  time  of  flight  for  a  parabolic  trajectory  may  be  derived  directly  from 

o  6 

the  equation  of  a  parabola  in  polar  coordinates,  namely,  r  =  r0  sec"  — ,  obtained  by  setting  e  =  1 

y  -|_  g  J 

in  the  conic  equation  r  =  — 2 -  (Figure  1). 

1  +  e  cos  6 


y 


Figure  1 .  Parabola  whose  focus  is  at  the  origin  and 
whose  vertex  is  at  O0,0) 
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We  have 


lfp 


J  6  hi  hi 


2  ,  -  rn  r  4  o  ,  ^  rft“  f  ->  0  ^ 


IT  =  I  r~dd  =  —  Isec4  —  00  =  —  Isec2  — 


h 


f 


0^ 


1  +  tan"  — 

V  2y 


00 


f  e2  e2  3 
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h  =  r0v0  =  r(u  ^  =  ^2/J^ 


2jur.  n  .  0,,  1  25' 

—  1  2  sin—  =>f,=-  — 2 

So  2  *  3V  M 


I -  eR 

r\  +  r2  +  Vrir2  C0S 


which  is  Equation  (17). 


DETERMINATION  OF  TYPE  OF  TRAJECTORY 


Proposition:  If  r,  and  r2  are  given,  then  tfh  <  tfp  <  lje 
Proof:  First  consider  Equation  (4),  namely. 


ffe  ~ 


(q  +  r2  )(.r  esc  x  -  cos  x) 


It  is  clear  that 


i\  +  r2  -  2 yfrfc  cos  v cos >  r{  +  r2  -  2  Jr\r\  cos ~^>  ri  +  r2~  2^[r\r\  cosh  y  cos ^ 


0, 


0, 


s  >  s0  >  s'. 


(18) 
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3 

d_ 

dx 


esc2  x(l  -  xcot  x)  ■ 


jcsc’  x(xcscx-cosx)}  =  csc;c{3csc2  v(l  -  .rcot  x)  - 1}>  0,  0  <x<x . 


(19) 


Hence,  esc2  x(.rcscv-cosv)  increases  monotonically  in  the  interval  [0,7r).  Since 

2 

limjcsc2  x(xcsc  x  -  cos  .r)|  =  — ,  it  follows  that 


esc2  v(.vcsc.v-cosv)  >  — ,  0<x<tt. 
v  ’  3 


(20) 


Observing  Equations  (18),  (19),  and  (20)  in  conjunction  with  Equation  (4),  we  conclude  that 
tje  >  t  fp .  Now  consider  the  expression  for  l  fh  derived  on  page  15,  namely, 


S' 


^  _  Q 

— — csch2y  <j  (r,  +  r2)(cosh  y—y  csch  y)  +  2^r\r2cos^-(y  coth  y  - 1)  f> .  (21) 


fh 
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We  have 
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=  — csch4y 
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3^2y  +  ^-  +  ... 
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3  1 


■2yU  +  ^r-  +  ..l-4y 
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n= 2 


(2/7  +  1)!  (2/7)! 


i2  y) 


2n+l 


=  -csch4yjT- 


:(2y)2n+3  <0. 


(2/7  +  3)! 

Hence,  csch2y(y  coth  y  -l)  decreases  monotonically  in  the  interval  [o,°°).  Since 

l 
3 


lim{csch2y  (y  coth  y  — 1)}  =  — ,  it  follows  that 

y— >0  ^  J  ^ 
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d_ 

dy 


csch2y  (ycoth  y-1)  <  — ,  y  >  0. 

jcscfry  (cosh  y  -  y  csch  y )}  =  csch  y |3csch2y  (y  coth  y  - 1)  - 1}  <  0. 


(22) 


Hence,  csch’y  (cosh  y  -  y  csch  y)  decreases  monotonically  in  the  interval  [0,°o).  Since 

2 

lim{csch2y  (cosh  y  -  y  csch  y)}  =  — ,  it  follows  that 

y-»0  1  '  "  J  3 


csch2y  (cosh  y  -  y  csch  y)  <  — ,  y  >  0. 


(23) 


Observing  Equations  (18),  (22),  and  (23)  in  conjunction  with  Equation  (21),  we  conclude  that 
tfh<tfp •  Hence, 

ffh  <  ffp  <  ffe  ■ 

Therefore,  given  rp  r2,  and  tf,  one  can  immediately  determine  the  type  of  trajectory  by 
computing  tjp  from  Equation  (17)  and  comparing  it  with  t  f  . 
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APPENDIX  A 

TIME  OF  FLIGHT  BETWEEN  TWO  POINTS  OF  AN 
ELLIPTICAL  TRAJECTORY 


A-l/A-2 
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For  an  inverse- square  gravitational  field,  given  r,,  vc,  and  |r2|  at  arbitrary  points  P]  and  , 
it  is  possible  to  derive  a  closed-form  expression  for  the  time  of  flight  between  Px  and  P2 ,  namely, 


where 


tfe=^-{E2~E  1  -  4sin  E2  -  sin  Ex )}, 
h 


-l 


h 


1  1  h 

(  2 

vn 

rixM,  P  =  ~, 

a  = 

— 

M) 

1-— ,  b  =  Jap, 


£j  =  arc  tan  VfN,  E2  =2  n  -  arccos 


1-* 

V  a) 


h  ( a  - 1] ) 

It  is  a  well-known  fact  that  the  focus  of  the  elliptical  trajectory  is  the  center  of  the  earth,  and 
this  can  be  taken  as  the  origin  of  the  rO  coordinate  system  (Figure  A-l). 


Figure  A-1 .  Position  and  velocity  of  an  object 

undergoing  elliptical  motion  with  the 
center  of  the  earth  as  the  focus  of 
the  ellipse 

The  derivation  is  facilitated  by  introducing  the  so-called  eccentric  anomaly,  E.  Choose  the 
center  of  the  ellipse  to  be  at  the  origin  of  the  xy  coordinate  system,  and  draw  the  circle 
circumscribing  the  ellipse  (Figure  A-2).  The  radius  of  this  circle  is  clearly  a,  where  a  is  the 
semimajor  axis  of  the  ellipse.  If  r  is  the  magnitude  of  the  radius  vector  at  P  and  6  the  angle 
between  r  and  the  x  axis,  then  r 2  =  (x-aeY  +  y2,  where  e  is  the  eccentricity  of  the  ellipse.  If 

the  perpendicular  from  P  to  the  x  axis  cuts  the  circle  in  (x,  y') ,  then  since  the  equation  of  the 

2  2 
X  V 

ellipse  is  —  +  -E-  =  1 ,  where  b  is  the  semiminor  axis  of  the  ellipse,  and  the  equation  of  the 

a  b 

2/22  y  y  y  b 

circle  is  x  +  y  =  a  ,  we  have,  for  the  same  value  of  x,  1  -  — y  +  —  =  1  =^>  ^  =  — . 

a  b  y  a 
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Differentiating  the  expression  for  cos  0  with  respect  to  time,  we  obtain 
-0  sin  0  ■ 


sin  E 1 

(l-  CCOSi?) 

l  -  eE  sin£’(cos£’-e) 

( 

is) 

1 E  sin  E  b2 
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[l  —  c  cos  E  ] 
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(1-CCOS.E) 

2  2 

1  r 

b  •  b-  ■  rO 

=>  —  0sinE  =  —  EsinE  =>  E  =  — . 
r  r”  Zz 

Now,  conservation  of  angular  momentum  implies  that  r  0  =  constant  =  /z . 

•  h  h 
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ab 
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/e  z. 
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where  p  =  a(l-e2). 
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aeh  eh 

Since  r  =  rr ,  v  =  rr  +  r0().  it  follows  that  r  •  v  =  rr .  Hence, 


.  br-x 

E  =  arcsin - =  arccos  t 
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V  a) 


=  arctan 
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The  energy  equation  is 


where  £  is  the  total  energy. 


1  o  Llm 

—  mv - =  £  , 
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Now 


■  h  .  -dr  h  dr  ,2  2a2  h2 
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Putting  r  =  — ,  we  have 
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Thus,  the  energy  equation  becomes 
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Putting  0()  =  --,we  obtain  the  polar  equation  of  the  orbit,  namely, 


h2 


, .  ,  2sh  . . 

r/<M+,  1h — ~ — cos# 

'  \  /LTm 


Comparing  the  above  equation  with  the  equation  for  r  on  page  A-5,  namely, 
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we  have 
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This  completes  our  derivation.  One  can  determine  whether  the  orbiting  object  is  before  or 
past  apogee  by  testing  whether  r  ■  v  is  positive  or  negative.  If  it  is  positive,  the  object  is  before 
apogee,  and  if  it  is  negative,  the  object  is  past  apogee.  To  prove  this,  we  have  r  ■  v  =  rr .  Since 
r  is  increasing  in  the  region  0  <  9  <n  and  decreasing  in  the  region  n  <9  <  2k,  it  follows  that 
r  is  positive  in  the  first  region  and  negative  in  the  second  region,  which  proves  the  proposition. 


As  a  relevant  mathematical  digression,  the  author  derived  the  condition  that  a  circle  of 
radius  R  intersects  an  ellipse  whose  focus  is  the  center  of  the  circle,  as  in  Figure  A-l.  The 

,  R}  L  Ir  .  ,  ,  . 

1 - ,  e  =  ,  1 - ^  ,  in  which  case  the  overlapping 

a  J  V  a 


required  condition  is  |//|  <  1,  where  Tj  =  — 


area  is 


R2  arccos  y — —  +  ab  |arccos  tj  -  e^jl-Tj2 ) . 
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APPENDIX  B 

AN  ALTERNATIVE  DERIVATION  OF  THE  TIME  OF  FLIGHT 
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First,  the  following  compact  relation  will  be  derived: 
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where  0 D  =  arccos  (r,  •  r) ,  6 R  =  arccos  (r,  •  r2) .  From  Appendix  A, 
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which  yields  Equation  (B-l).  Now, 
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It  can  be  shown  after  some  algebra  that 
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This  is  Equation  (15)  of  the  main  text. 
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APPENDIX  C 

EXPRESSIONS  FOR  /,  g,  f,  AND  g  IN  TERMS  OF  rp  r2,  AND  p 
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Proposition: 
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where  0R  =  arccosft  ■  r2 ) . 
Proof:  From  Appendix  A, 
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-cos  — - 
2 


-  sin  — . 
P  2 


cos^- 
P  2 


rMl~e2)  .  0,-0,  4^2  .  0R 

-sin  — - -  =  — - sin  — 

b  2 


P  2 

sin  (2.r)  =  sin  (E2  -  E, )  =  sin  E ,  cos  E,  -  sin  E,  cos  E 


(C-l) 


=  -^-sin#2  — (<?  -i-cos#,)  -  — sin#,  -^-(e  +  cos#,)  =  — {e(sin#2  -  sin#,)  +  sin(#,  -0,)} 

bp  bp  bp 

2  r.r2  .  0,-0, 

-  — — sin^ - L 


bp 


0,  +  0,  0,  —  0, 

ecos— - —  +  cos  — - L 


2  r  v  f) 


bp 


■R  0,  +  0, 
cos  —  +  ecos— - - 


Or 

2 


Now 


rx+r2  = 


-  +  - 


=  P- 


2  +  c(cos#,  +  cos#9 )  2  r,r,(.  0R  0,+0,^ 

v  1  22-12  1  +  C  COS  — —  COS  — - - 

l  2  2  y 


1  +  ecos#,  l  +  <?cos#2  (l  +  ecos^Xl  +  ccos#,)  p 

0,  +  0,  0R 


ecos- 


f  p 

fi  0 

J 

-1 

l2 

lrl  r2j 

J 

Substituting  this  into  the  formula  for  sin(2.r) ,  we  obtain 

f  i  i  A 


sin(2.r)  =  sin  — 

Pp  2 


1  1 

—  +  - 

vri  r2  y 


■1 


cos 


0C 


0R 

-  +  cos  — — 


2r.r,  0K  p 
-  — —  tan  —  \  — 


bp 


2  2 


1  1 

— +  — 

Vri  V2J 


■  1  +  cos2  — 
2 


r,r,  0„ 

—  tan  — 
b  2 


1  1  2  .  ,0R 
—  + - sin-  — 


*i  r2  P 


V  'l 


r,r,  0K 

—  tan  — 
b  2 


fl  1  1  —  cos  0R  ^ 
- h - - 


Vri  r2 


(C-2) 


J 
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It  has  already  been  shown  (pages  4  and  5  of  the  main  text)  that 

QR  ■  Jjua  .  .  s  .  2a  .  2 
;  — =  — - - sin(2jt),  g  =  1 - sin  x. 


la  .  2 


/  =  1 - sin"  x,  g  =  2„ 


m-,r2 


sin  x  cos - 
//  2 


rir2 


Using  Equations  (C-l)  and  (C-2),  we  obtain 


2ar2  .  ,  6 

f  =1 - T^-srn 

b~  2 


2  0;„2  UR  _  y_.2  „;„2 


2  A  .  2 

-sin 

P  2 


g  =  2  I—  —  sin— cos—  -  , — 
>  *  2  2  ^ 


sin6> 


•  Jua  0K  f  1  1  1  -  cos  6,  x 

/  =  —  -^2- —  tan  —  |  — +  — 


Vri  E 


M  &R 

=  |—  tan  — 
2 


1  -  cos  1  1 

p  n  ri  j 


.  2 or.  .  n  0R  2r, 

y  =1 - J_sm--5-  =  l - i 


sin2  — 
P  2 


which  are  the  desired  results. 
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APPENDIX  D 

TWO  DIFFERENT  SETS  OF  INITIAL  CONDITIONS 
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Set  I.  Given  r, ,  r2 ,  and  vc  ,  determine  the  type  of  trajectory  and  its  solution. 
We  have 

/i  M  j 


.2  'N 


If  a  >  0,  the  trajectory  is  elliptical;  if  a  <  0,  it  is  hyperbolic;  and  if  a  =  °°,  it  is  parabolic. 
In  the  elliptical  case,  we  have 

S  =  rx  +  r2  -2^/rp^cosvcos^-  =  2a  sin2  v  =  2 a  (l  -cos2  v) 

Q 

=>  2 a  cos2  x  -  cos  x cos  +  rx  +  r2  -  2 a  =  0 


v  =  arccos 


—  <  cos  — ± .  r.n  cos2  —  -2 air,  +  r,  -  2a) 
2  a  s  2  V  2  v  ’ 


where  0R  =  arccos  (r,  •  r2) .  Thus,  the  elliptical  case  is  ambiguous. 
In  the  hyperbolic  case,  since  a  <  0,  it  follows  that 


0, 


2 a  (rj  +  r2  -  2a)  <  0  =>  J rxr2  cos2  —  -  2 a  (rx  +  r2-2a)  >  ^rxr2 


cos- 


6„ 


Hence, 


—  <  x[rr\  cos  —  -,  r.r^cos2  —  -2a(t:  +  r0  -2a)  > > 0, 
2a  V  1  2  2  V  2  V  1  2  M 


—  <  JjpTc os  — +.  r.r^cos2  —-2a(r,  +  r, -2a)  l<0. 
2  a  v  1  2  2  V  2  V  1  2  M 


Putting  a 


1 


r,r2  cos  — - . /r,r2 cos2  — -2a(r,  +  r2 -2a)  1  and  noting  that  coshy>l,  we 


2a  2  V  2 

surmise  that  w  >  1.  To  prove  that  this  is  indeed  the  case,  we  have 

8, 


rx  +  r2  >  2y]  rxr2  cos  =>  -2a  ( rl+r2)>  -4  a^t\r2  cos  - 


t  6  f 

rxr2  cos2  - 2a ( t\  +  r2  -2 a)> 


V 


eR  . 

rr,  cos  — -2a 
2 


Y 

) 


t\r2  cos^  ~^~2 a(rl  +  r2  -2a)  >  Jr\r2  cos -  2 a 


Q  _  /  , - 

rxr2  cos2  —  -2a(r1  +  r2  -2a)  -^)p^cos  —  >  -2a  =>w>l  =>  y  =  arccosh  w  =  ln|w  +  yju2  -1 


e 
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Set  II.  Given  r,,  vc,  and  tf,  determine  the  type  of  trajectory  and  its  solution. 
As  in  I, 


a  = 


2  v 


2\ 


-1 


_ _ C_ 


and  the  type  of  trajectory  is  determined  by  the  sign  of  a  .  Now,  from  Appendix  A, 

Z7  1  'i  ■  r  brr  V  r,  •  vc 
ecos Ex  =1-  —  ,  esin Ex  = — 1 — -=  ' — £ . 

a  ah  y Jjua 


Hence, 


=  a  (l  -  e cos  E2 )  =  a  {l  -  e cos  (2.x  +  £j ) }  =  a {l  -  e cos  El  cos  (2x) +  e sin  £j  sin  (2.x )} 


=  a  1  - 


f  r\ 

1  — — 


r,  -v 


cos(2.r)  +  1|/— sin(2.r)  J >  =  a-(a-  rI)cos(2x)  +  1—^  ■  vc  sin(2.r). 


V  aj  JJm 

From  Equation  (1)  of  the  main  text, 


A 


tfe  =  2,1—  {x  -  e  sin  x  cos  ( x  +  £j )}  =  2  1—  ( x  -  e  sin  x  cos  x cos  £j  +  e  sin2  x  sin  £j ) 


M 


M 


=  J—  |2x-  e sin  (2x)  cos  Ex  +  2e,sin2  xsinEj} 


\a-\2x- 


(  r\ 
1--4- 

l  a) 

dt 


.  ,  v  2rrve  .  2  , 
sin(2.r)H - |^=^sm  x 


J/ua 


—  {2 ax  -  (a  -  rj )  sin(2x)}+  —  Tj  •  ve  sin2 


la 


M 


dx 


~  =  2J-  {a  -  (a  -  ri  )c°s  (2x)}  +  ^-r,  ■  vc  sin  (2x)  =  2 x\—r2 . 


M 


I  a 

M 


Referring  to  pages  4  and  5  of  the  main  text,  r2=  frx  +  gvc,  r2  =  /r,  +  gvc ,  where 


/  =  1  -  —  sin2  x, 

ri 

g  =  ^{sin  iE2  ~  E\ )  -  e  (sin  E2  -  sin  Ex )}  =  tfe 


-{2.r-sin(2.r)}, 


•  yj /.la  .  \  ,  2(3  .  2 

/  = — - - sin(2x),g  =  l - sin  x, 


where  x  is  updated  by  the  Newton-Raphson  method  as  follows: 
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ltf~tfe(Xn)_  ,  1 

Xn+l  ^  '  X" 


Jt 

v  dx  J 

\  /  x=xn 


2riM 


tf-tftM}’ 


and  a  good  initial  guess  is  x,  =  0. 

The  corresponding  equations  for  a  hyperbolic  trajectory  are: 


r2  =  a-(a-rl)cosh(2y)  +  r \  ■  vcsinh(2y), 

tfh=^~f  {2ay  “  (a  “  ri  )  Sinh  (2^  )}  ~^'fc  Sinh2  *  ~^  =  1\ 

/  =l  +  y-sinh2  y,  g=tfh  +  aj^-{sinh(2y)-2y }, 
f  =-^~^a-sinh(2y),  g  =1  +— sinh2  y, 


where  y  is  updated  as  follows: 


y«+i  -  —  -  y» + 


fh 


V  5>’  J 


y=y„ 


and  a  good  initial  guess  is  y,  =  0. 
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APPENDIX  E 

DERIVATION  AND  COMPUTATION  OF  THE  ELEMENTS 
ARISING  IN  THE  0-GUIDANCE  MATRIX 
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Suppose  that  the  position  vector  at  release  is  slightly  perturbed.  Let  r0  and  r,  denote  the 
unperturbed  and  perturbed  position  vectors  at  release,  respectively,  and  r2  the  position  vector  at 
impact.  Then 

9ro  =  arccos  (?0  ■  r2 ) ,  0O  =  r2  esc  0RO  -  r0  cot  0RO,  n0  =  r0  x  0O,  r,  =  -  (rlrr0  +  rig%  +  rlp n0 ) , 


r,  =  ,/n2.  +  r~e  +  r2  ,  dR  =  arccos  (r,  ■  r2) ,  =  r2  esc  0R  -  f  cot  0R. 


_  rlr ,  =  _riP 

drlr  /-  '  3rw  rt  ’  drlp  r\ 

dr,  1 


A  a-* 


=  1, 


= r'2)po  - (riA + vf,)} 


drlr  i] 


—  =  4{(ri2  -  4)4  -  rw  4>o  +  riTo)} 


dr, 


\e  'i 

dr,  1 


af  -  =  J  {('i2  -  ■ r‘p)  fl»  “  V  +  rj„ )} 


dr,P  n 


f  dr,  ^ 

f  df, A 


=  0. 


=  0, 

Gn 


v3rw^0 

^  dr,  ^ 


f  A  /v 

4f_  =£ii 

,  dr.  rn 

v  i/>  yri=r0  ° 


Differentiating  cos  0R  =  r,  ■  r2 ,  we  have 


dr. 


R 


■  sin 0  ^-  =  -^4f 

5111 UR  -v  *2 

drb.  drlr 

CSC0B 


dr. 


lr 


rl 


f/  2  2  ' 

i  * 

£  ( 

1 

1 

If. 

'  L  -  r\r  ( 

CSC0B 


{(4  -  4 ) cos 0RO  -  rh.rw  sin  0RO) 


=  0; 


fd0,  2 


30, 


R  _ 


dr, 


=  —csc0r 


dr, 


id 


dr. 


L  = 


CSC£L 


{(4  -  4  )0(,  ■  4  -  rw  (ri>0  •  4  +  rj0  ■  4 )} 


10 


esc 


{(4  -  4 ) sin  4o  -  A  A  cos  dA>0} 


4^4 


d6>, 


/?  _ 


dr 


=  -esc  9. 


df, 


ip 


dr. 


r,  = 


CSC  6L 


ip 


{(4  -  4 )  4  ■  4  -  'ip  ('iTo  ■  4  +  r„0()  •  4 )} 


r{  esc  9 R  .  . 

=  - (A  cos  9ro  +  rw  sin  0RO ) 


fd_9^ 


=  0. 
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=  ~r2  esc  0R  cot  dR 
dr.  dr. 


-  +  fj  CSC2  0R  <^-!L  _  -^-cot  0R 


1  r 


dru.  3rlr 


dO  3r 

:  csc#R  (r,  csc#fi  -r2cot#fi)  — --cot#  1 


dru.  ~~“'*3rlr 


02  esc  0R  — -  -  cot  0R  -^-L  =  esc  0R  (r,  sin^  -0,  cos  0R )  ^  -  cot  6 

nr  nr  '  '  nr 


3r, 

drh.  K  drh 

=  (fj  -0,  cot  0R  |  — — ^ —  cot  0R  -^L 
V  1  1  R)  R 


3  0„ 


fd^ 


=  0. 


dr. 


dr, .. 


Similarly, 


30 


^  =  (f1-01cot^)|^-cot^j 
orle  orw 


dr, 


3tl 

dr. 


\e 


(r0  0ocot0so) 


cotggoA  __r± 

“n  ’ 


ro 


ro 


30, 


[j  /  A  \  ()  0 

—  -  (fj  —  0J  cot  ^ —  cot  0 

i  p  p 


dr. 


dr. 


dr, 


i  p 


cot  0 


RO 


ro 


no- 


dvCj  _  dvc,  dr,  |  dvci  d0R 


drh.  dr,  dru.  d  0R  d rlr 


dva  _  dva  ch\  |  dvc,  d0R 
drw  dr,  drw  d0R  drw 

dvci  _  dvci  dr,  3vd  d0R 


A=r, 


^3v  A 


rdv^ 


y  dr,  _ 

-r0  V  i  /rt=r0 


rdv  A 


v  d0R  )_  _ 
V  R  Jr,= r0 


r  dv  .  A 


dr,p  dr,  dr,p  d0R  dr,p 
Differentiating  vc  =  vcrr,  +  vcj90, ,  we  have 


vsVy„ 


=  0. 


3vr  3vrr „  3r,  dvr0  -  30, 

a — — r'  +Vcr^~  +  li — 01  +Vcdz~ 
3  r„.  3  r1(.  drh.  3  rlr  drh. 

3v„  3v 


r 

f  3v  ^ 

f0  + 

r,=rn 

l  3ri  j 

f,=rn 

V  dri  J 

en; 


3r.  3r  r‘+Vcr 


5?1  ,  ^,0  A  ,  ,,  39| 


^3v  ^ 


V3rWfl=r0 


'10 


10 

'3v_  A 


"  +  "^7^0i  +vce  jr 
or,  q  or, 


r\G 


\d0Rj_  _ 
V  K  /ii=*b 


+  V 


cB  0 


1 


lr0 - \ 


10 


^3v  A 


cB 


\d6Rj-^ft 


-  V 


crO 


0 


0’ 


3vr  3vrr „ 

VJL  =  T^ri+v«- 
dr,  „  dr, 


'i  p 


'i  p 


3f,  ,  3v,„a  ^  3^ 
3r; 


^  +  T^0l+Vf0 

dr,  „  dr, 


'i  p 


'ip 


ip 


rd\  A 

v3Va=,, 


(VerO  Vc0O  C°t  )  A0  . 
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It  is  desired  to  evaluate  the  partial  derivatives 

Differentiating  the  formulas 

2 jur,  .  0R  \2ju  (  \*>  &R 

VcS=t— Vcr  =  \—  J^COS^f-C 

2  V  S 


^.1  ( and  f dv" 

drl  drl  J,/  ld6>R 


we  have 


6 R  \  )2{X  (  ITT  9 R 

— -cosx  ,  v,  =.  — —  cos.x-  —cos  — 
2  "VS  Vr.  2 


r,  sin  x  sin  9,, 

hs3  - 


I  sin  —  [  2r\  +  r2  -  3  Jr,  r2  < 
1 2 r3S3  2(  1  2  V  1  - 


r,  r,  cos  xcos — 
-  2 


I  Llr S  f/  \  dr,  I -  9  9„ 

=  - =-r1  (h  +  ^  Icos  — -,/r,r,  cosx  1  +  cos'  — 

y 2rj53  [V 1  2  2  Vl2  ^  2 

/  2//  f  ■  2  @r  I -  @r  ' 

-  =  J— Vsinv  r, +r0sin  — 1 -  -  ,/r,r,  cosx  cos—1 -  , 

V  S3  l  2  *  2  J 

//£'  .  tVy,  f  .  9  /  tV/, 

=  -  ~  sin  — -  r,  sin'  x  +  r0  -  Jr, r0  cos x cos  — 

A  2r,53  2  l  1  -  V  i  -  2 


2//  .  f  ■  2  @R  I - 

— V  sin  x  r,  sin  —  +  r,  -  Jr, r,  cos  x cos  — 
S3  ^  2  -  v  i  2  2 


I  i  r  *  2 

*  sin  — -  r,  +  r?  sin  x  - 
j  2r2S3  2  l1 


(E-l) 

(E-2) 

(E-3) 

(E-4) 

(E-5) 

(E-6) 
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Now 


+  dvce 

(  dx  " 

( 

(  dvc^ 

+  dvcg 

f  dx  ^ 

1 6rl  J 

1 6rl  J 

dx 

X 

l  drO 

lf 

l  Mr) 

lf 

l  Mrj 

dx 

X 

with  similar  expressions  for 
derivatives,  we  have 


and  I'3’'"  '' 


.  Using  a  well-known  theorem  in  partial 


dx 

\drJ, 


dr, 

A  j 


d_!i 

v  dx  J 


-1 


rdx} 

v5riy 


V5ri  J 
f  dtj,'] 

Kdxj 


'v 

V3ri  J 
dtf 

dx 


dtf 

where  — —  is  given  by  Equation  (5)  of  the  main  text  for  elliptical  trajectories,  and 
dx 


fc ^ 
v  a'i  y 


CSC2  X 


lyjl/Is 


Q 

3  (t\  +  r2)  (xcsc  x  -  cos  x)  -  6  r2  cos2  cos  x  (l  -  xcot  x) 

W  o 

+  —  cos— {3(3^  +  r2)(l -xcotx)  ~{5t\  +  r2)sin2  xj 
V  r,  2  L  '  J 


obtained  by  differentiating  Equation  (4)  of  the  main  text.  Similarly, 

fdt,  ^ 

r  dx  A 


\d0R  j 


\d&R  j 


dt 


f 

dx 


^dt  ^ 


where 


v  d  6  R  j 


is  given  by  Equation  (6)  of  the  main  text. 


Using  the  above  results,  It  is  straightforward  though  tedious  to  verify  the  following  three 
identities: 


rdv„  A 


+  vce  +  rx 


V  dri  J 


=  0, 


(E-8) 


r  dv„  A 


\  dOR  j 


+  V 


c9 


\d&R  j 


+  r, 


dv 


c9 


1  dt 


0, 


fdv  ) 

cr 

r  ^tr  ^ 

\^9r  j 

ri 

k^9r  j 

f 

3  dv 


+  ^ce=~tf 


C0 


2  dt 


f 


(E-9) 


(E-10) 
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From  Equation  (E-8),  we  have 

dv 


cG 

dr{ 


Jt' 


f  3..  A 


dv 


\d0R  jt 


+  v 


cG 


(E-8') 


which  proves  the  symmetry  of  the  (^-Guidance  matrix. 

3v„,  ' 


Eliminating 


\^0R  Jt 


between  Equations  (E-8)  and  (E-9), 


dv 


cG 


\  drl  Jt 


cG 


dv 


cG 


Jt 


-  V. 


+■ 


3v 


cG 


dt. 


(E- 11) 


Writing  Equation  (E-8)  for  the  reverse  trajectory,  we  have 


f  A 


3v 


\^°RJt 


( A 


to 


dv 


V  dr2  J, 


0. 


Now 


r\ 

vw  =  ~vce’ 


( 3..  A 


to 


dv 


V  y, 


cO 


dv 


V  y, 


c0 


(E-12) 


A  A 


3v 


V^OrJ, 


fdv  A 


:  E*  +  E 


v  f>2  yf 


A  :v.  A 


3v 


v  y, 


(E-13) 


Using  Equations  (E-8)  and  (E-13)  to  eliminate 


^3v„  A 


J 


and 


f  dv..  A 


\d°Rj, 


from  Equation  (E-10), 


^3vc/ 


V  dr2  Jt 


2  V°e  +  Vx 


f^vce^ 


3q 


3  3v 


+  -t 


cG 


f 


Jt 


dt 


f 


(E-14) 


Writing  Equation  (E-l  1)  for  the  reverse  trajectory, 


V  *2  j, 


_  VtG 


\d0R  Jt 


+  V , 


+  ■ 


dv 


tG 


dt. 


Substituting  the  system  of  equations  (E-12)  into  the  above  equation, 


3v 


cG 


V  dr2  Jt 


r0v 


2  vtr 


v  cG 


dv 


cG 


\^0r  Jt 


+  r. 


dv 


cG 


dt. 


(E-15) 
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Eliminating 


V  dr2  J 


between  Equations  (E-14)  and  (E-15), 


\ 


c8 


r,v 


dv 


d  0L, 


2vtr  \"WR  J 


2Vc0  + 


^  r2  3  dv 
— - 

VE,  2  j 


c8 


dt 


a 


c8 


f 


dv 


V  dri  J 


=  0. 


Eliminating 


f  3..  A 


c8 


dv 


V  ^  J 


between  Equation  (E-ll)  and  the  above  equation, 


V  ) 


^crVce-y ,  -Vcr 


To  3 

— - tf 

\vtr  2  y 


dv 


c9 


dt , 


c9 


f  r,v  A 

_  '  1  v  cr 

v 


(E-16) 


The  expression  (E-16)  is  indeterminate  when 

rtvcr  =  r2vtr  =>fr\c=f2-\t  =^>  sin  Ei  =  sin E2  =^£,  +  £,  =  3^ 


(see  Appendix  A),  which  can  occur  when  El  is  in  the  third  quadrant  and  £2  is  in  the  fourth 

quadrant.  Therefore,  it  is  desirable  to  derive  an  expression  for  the  above  partial  which  is  analytic 
everywhere.  This  is  done  in  Appendix  F. 

The  remaining  partials  in  the  Q  -Guidance  matrix  are  now  readily  obtained.  Employing  the 
equations  for  vcr  and  vtr  in  terms  of  vc9 ,  namely, 


we  deduce 


v„„  = 


cot#s  -  —esc  9r 


Vc8  + 


-tan- 


J 


riVc8 


V,r  = 


A 


esc  6r  -—cot  0R 


J 


A  0R 

vcd - tan^  = 

r,v~  2 


f  r  A 

1-^ 

V  hj 


0„ 


VC8  cot-f- v  , 


(dv„  A 


\^0rJ 


-VlrCSC0R  + 


-A - V 

1  < 

K 

1 

(  r \ 

1  +  -3- 

tan—} 

l  r2j 

2J 

,  ~dOR  , 

(E-17) 


V  drl  J 


A..  A 


dv 


Jt 


+  v 


c9 


(E-80 
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foe, 

dr{ 


J, 


f  r  ^ 

1  +  -5- 

V  r2j 


Q 

tan  —  +  cot  S-v 
2 


/? 


+ 


V,r 

VC9 


f  r  ^ 

1  +  i 

v  r2y 


On 

tan  — 


dv 


c8 


dr. 


Vtr 

VC0 


l  +  i 

V  r2j 


0R 

tan  —  >< 


vcd 


+ 


dv 


C0 


dr{ 


h. 


' — vc6  cot  eR 


[r  vf 

(  rA 

eS\ 

(dv  \ 

1 

1+^ 

tan  —  > 

cr 

[vce 

l  r2  J 

2J 

A®rJ, 

+  vce  c°t  QR 


where  Equation  (E-8)  has  been  used.  Substituting  Equation  (E-8'j  into  Equation  (E-14), 


dv 


cB 


V  dr2  Jt 


2  Vc0  + 


dVg 

J, 


3 

—  t. 


dv 


cB 


dt , 


Finally, 


<K 

V  dr2  Jt 


\vce  csc^+p— 

h  \VcB 


1+* 

V  r2j 


e R 

tan  — 


dv 


cB 


\  dr2  Jt 


E-9 


(E-18) 


(E-19) 


(E-20) 
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APPENDIX  F 


DERIVATION  OF  AN  ANALYTIC  EXPRESSION  FOR 


Vi 

)t 


F-l/F-2 


NSWCDD/TR-03/8 1 


As  stated  in  Appendix  E,  the  expression  (E-16)  is  indeterminate  when  f\vcr  =  r2vtr.  An 
( 3v  ^ 

analytic  expression  for  — —  will  now  be  derived  in  two  ways. 


First  method:  From  Appendix  A, 


rx  =  tf(l-ecos£’l),  r2  =  a  (l  -  ecosE2) . 


From  Equation  (C-l), 


0R  b 

sin  —  =  ,  sin  x  . 


From  page  2  of  the  main  text, 


0„  a 


{cos  x  -  ecos  (x  +  Ex )}. 


0„  6„  lab 


_  .  Up  Up  Z.UU  r  /  \1 

0R  -  2 sin— cos  —  = - smxjcosx-  ecos(x  + 


2  2  rxr2 


-{sin  (E2  -  Ex )  -  e  (sin  E2  -  sin  Ex )} . 


cos  6r  =1-2  sin2  —  =  1  -  sin2  x 
2  rxr2 


-{(1  -ecos  Ex )  (1  -  ecos  E2 )  -  2  (l  -  e2 )  sin2  xj 

-{cos  (E2  -Ex)-e  (cos  Ex  +  cos  E2 )  +  e2  (1  -  sin  Ex  sin  E2 )} . 


_  0  W  _  0 

K  +  ri~  2riri cos  6r=  rx+r2  +  2 7^ cos  -f-  rx  +  r2  -  2^  cos 


=  a 2  {2  -  e  (cos  Ex  +  cos  E2 )  +  2 cos  x  -  2ecos  (x  +  Ex )}. 
{2  -  e  (cos  Ex  +  cos  E2 )  -  2cos  x  +  2ecos  (x  +  Ex )} 

=  4a2  sin2  x{l  -  e2  cos2  (x  +  Ex ) J . 


From  Appendix  A, 


_rn  m 

vc6  ~  '\u\  ~  ~  ~  \ 
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.  ah  .  JJui  . 

=  rx  =  — esinEj  =  — - esinE, , 

(F-8) 

br{  rx 

vlr  =  r2  =  ^~  esmE2. 

(F-9) 

Substituting  Equations  (F-3)  and  (F-7)  into  Equation  (1 1)  of  the  main  text  and  using  Equation  (F-4), 

\I~F  =2  j]  ^Jr^vle  cos  —  sin  x  =  sin  x  {cos  x  -  e  cos  ( x  +  El ) } 

2  r, 


=  ^-{sin  (E2  -  E, )  -e(sinE2  -sinE,)} 


sin  0R  =  —  Ve  . 
jubr2 


(F-10) 

(F-4') 


From  Equations  (F-2),  (F-6),  (F-7),  and  (F-10) 
(r,2  +  r2  -  2rjr2  cos  0R )  -  4 /rr22  sin 


2  4 
ri  vce 


Or  4//  Z?  2  „r _ 2  ..  „2 _ 2 


sin  xicos"  x  -e  cos 


(x  +  Ej} 


L: 


|sin2  (E2  -  Ej )  -  e2  (sin  E2  -  sin  E1 )  | 

#1 

—  {sin  (E2  -  Ej )  +  e  (sin  E2  -  sin  E, )} . 


(F- 11) 


Also, 


tf  =  — {E2  -  Ej  -  e(sin  E2  -  sin  Ej )}, 


A 


(F- 12) 


arctan  ■ 


r\  fa  +  r2  )Vc8  ~ 


E1  —  E 

=  x  =  — - -  (Equation  (12)  of  the  main  text). 


Substituting  Equations  (F-2),  (F-4'),  (F-7),  (F-ll),  (F-12),  and  (12)  of  the  main  text  into 
Equation  (14)  of  the  main  text,  namely, 

dtf  1 


-  ~  — J  (ri  +  r2  -  2rir2  cos 0R )  -  4 /rr22  sin4 


dvc0  F 


h  +4^3E3v>in3^arctan_ 


Ve 


V 


eVe 


ri  (ri  +  0 )  vle  ~  2Fri  sin2 


\6fU2r2r2v2cg  sin3  0R  sin 


2  9 R 
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we  obtain 


dtf  a2b 


dvrff  Vf 


{sin(F2  -Fj)  +  e(sinF2  -sinF,)}. 

{3(F2  -F,)-e(sinF2  -  sin  F,)}- 16  sin2 


(F-13) 


Substituting  Equations  (F-7),  (F-8),  (F-9),  and  (F-12)  into  Equation  (E-16),  namely, 


1 


dv 


VcrVc9 


r ,  -  v„ 


c9 


j, 


\Vtr 


dv 

dt. 


c6 


vc9 


'  r,v  " 

^  _  1  cr 


V 


r2Vtr  ) 


we  obtain 


v  3^*  /, 


2r2vcS  (sin  F2  -  sin  F, ) 


/ube  sin  F,  sin  F2 

2  ( r,2  sin  F2  -  r22  sin  F, ) 

+3 a2e  sin  F,  sin  F2  [F2  —E1—e  (sin  F2  -  sin  F, )  j 


dv 


c8 


d  t, 


From  Equations  (F-l), 

r2  sin  F2  -  r2  sin  F,  =a2  j(l  -  <?cos  F, )~  sin  F2  -  (1  -  ecos  F2 )“  sin  F,  J 

=  a 2  (sin  F^  -  sin  F, )  |l  +  e2  ( 1  +  sin  F,  sin  F2 ) }  -  2e  sin  ( F2  -  F, ) 

=^>  2  (r,2  sin  Ft  -  r2  sin  F, )  +  3  a2e  sin  F,  sin  F2  {F^  —Ex—e  (sin  F.,  -  sin  F,  )J 
=  a 2  |jsinF2  -sinF,  ){2  +  e2  (2  -  sinF,  sin  F^)|+  e\3 (F2  -  F,)sinF1  sinF.,  -4sin(F2  -  F,)}J. 

d 


dv 


C0 


2/,2v((9  (sin  F2  -  sin  F, )  d/  ( 


//de  sin  F,  sin  F2 

dv 


/ 


rtf 


-a2r\  <; 


|^(sinF2  -sinF,  )[^2  +  <?^  (2-sinF1  sinF2)^| 
+e[3(F,  -F1)sinF1  sinF2  -4sin(F2  -F,)J 
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Substituting  Equations  (F-10) 


ana  (r- 


f  3v  ^ 


vc0 


d0R  i 


fua 

^r\vce  (sin  E2 


' 2b 2 


E2  -  sin  Ex 


_ <Kl 

,)Vf  dtf 


linE’,  -sin£’1)]-16sin2 


! uab 2  dvc0 
rivce^F  dtf 


"(sin  E2 


[sin(£’2  -Ei)  +  e( si 

esmExsh\E2\ 

[3(£’2  ~E{)  -e(sin 
-{sin(£’2  -  Ex )  -  e (sin  E2  -sinE,)}. 

^(sinE^  -  sinE,  )[2  +  e1  (2  -sinEj  sinE2)^| 

-  E, )  sin  Ex  sin  E2  -  4sin  (E.,  -  Ex  )J 


+e[3{E: 

-sinEj^  +  e2  -e2  sinEj  sinE2)-(3e2  +  l)sin(E2  -E,) 

+3e2  ( E2  -  E, ) sin E,  sin E. 


Using  Equation  (F-7)  in  the  above  equation, 


V  J 


dvc0  j e(sin  E2  -  sin  E{  )(3  +  e2  -e2  sin  Ex  sin  E2 ) 


(  4f  dtf  j-(3e2  +  l)sin(E2  ~El)  +  3e2{E2  -E^sinEj  sinE2 


•  (F-14) 


Note  that  a  and  b  can  be  computed  from  Equations  (3),  (9),  (1 1)  of  the  main  text,  and  (F-7).  We 
have 


a  2  Scscr  v  2  «  R,b  rlVJ 


r2r2v20  sin  0R 


Second  method:  Equation  (F-14)  may  also  be  derived  from  the  identity 

=  -l. 


{  dvc&  ] 

fdSK' 

(  ^1/  ) 

J 

h 

l3'/  J 

VcH 

Kdvc0j 

Or 

'dt  ^ 


f 


\dvcej. 


is  given  by  Equation  (F-13).  Differentiating  Equation  (13)  of  the  main  text,  namely, 


r,r2vrg  sin  9,. 

tf  ~  1  .  2 Hrxryce  sin"  9j 


Q 

2 Pr2  ( rx  +  r2 )  sin2  -  rxv2e  ( r~  +  r2-  2 rxr2  cos  9R ) 


+  - 


Ve 


arctan- 


Ve 


ri(ri  +  r2)vc0-2//r2sm2^ 


we  obtain 
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fV 

f 

—  t 

~  '/ 

v 

cot  6,, 


1  dF 


F  d0 


+ 


Rj 


rir2vce  sin2  dR 


V(rl  +  r2)-2ryc0  +  ^Jism2^ 
r,v„a  2 


+  - 


LI  0 

•  +  I 


+ 


2  rxvcgF  2 

lir\r2vle  sin  #A, 


tan ^-\rAri  +  ri)  vle  -2//r2sin2^~ 


4cot#„ 


1  dF 


V 


arctan  - 


s[F 


R  J 


rAn  +  r2)Ke-^Mr2^2Y 


(F-15) 


Differentiating  the  expression  for  F  on  page  1 1  of  the  main  text,  namely, 


0, 


F  =  \2Mr2  shr^-rpy2. 


r 


0 


r1  +  r2-2yJrlr2cos  — 
v  z 


riKe 


n  +  r2  +  2^2 cos  ■y  J _  2  ^2 sin2  y  |  ’ 


X 


,2  Qr 


we  have 


|y  =  2/2  sin  y  j-r/v^  +  //rt  (rt  +  r2 )  v2e  -  2fi\  sin2  M . 


Substituting  Equations  (F-l),  (F-2),  (F-4'j,  and  (F-7)  into  the  above  equation, 
-3y  _  2 flab  yy |cos^^  -  £) )  -  <?(cos  El  +cosE2)  +  <?2}. 


dOR  rt 


(F-l  6) 


From  Equations  (F-l),  (F-2),  and  (F-7), 

//(r,  +  r2)-2r^v2cg  +  ^-p-sin2  —  =  fua  {l  -  cos  ( E2  -  Ex)  -e(cosF,  +cosF2)  +  2<?2}.  (F-17) 

riVc0  2 

From  Equations  (F-2)  and  (F-4'), 


Or 

tan  —  = 


sin 


0R 
2  _ 


cos 


0* 


2 sin2—  _  ,3  •  2 
2  _  2/jb  sin  x 

sin  0R  ar{  4f 


(F-l  8) 


From  Equations  (F-7),  (F-16),  (F-18),  and  (12)  of  the  main  text, 

yjyy tan y  |ri  ( ri  +  h ) Ke  ~  2FF  sin2  y } |y  =  yj=T sin (2x) {cos (E2  -  El )  - e (cos Ex  +  cos E2)  +  e2} 
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Substituting  Equation  (F-10)  into  the  above  equation, 


eA  dF 


~TF tan  j  ri  ( ri  +  ri )  vle  ~  sin2  yj  })Q 


2  i\v-ceF  2 


jua  sin(E2  -  Ex) 


sin  (E2  -  El )  -  e  (sin  E2  -  sin  E, ) 
From  Equations  (F-4'),  (F-5),  and  (F-16), 


{cos(E2  -  E, )  -  e (cos E,  +  cos E2 )  +  e2 }. 


(F-19) 


4cos^  -  — ^—smOR  =  -^-{cos(E2  -Ej  -e(cosEt  +cos E2)  +  e°  (l  - 2sin  E,  sin  E2)}. 
F  oOr  rxr2 

Dividing  the  above  equation  by  sin  0R  and  using  Equation  (F-4), 


4  cot  0n 


1  dF 


2a- 


F  ddR  i\r2  sin  0K 


{cos  (E2  -Ex)~e  (cos  E,  +  cos  E2 )  +  e2  (l  -  2  sin  Ex  sin  E2 )} 


_  2 a  cos (E2  -  Ex )  - e (cos Ex  +  cosE2)  +  e2  (l - 2sinE,  sinE2) 
b  sin  (E2  -  E, )  -  e (sin  E \  -  sin  E, ) 

From  Equations  (F-4'),  (F-7),  (F-20),  and  (12)  of  the  main  text, 


(F-20) 


Fr\  r2vle  sin  f 


Ve 

Ma  ( E2  ~  E| ) 


4  cot#,-, 


1  dF 

~FdK 


arctan  - 


Ve 


eD 


ri  (ri  +  h )  Ke  ~  2//r2  sin2  -y- 


V  1  ) 

cos (E.,  -E,)-e(cos E,  +cosE2)  +  e2(l-2sin  E,  sinE.,) 


(F-21) 


sin  (E2  -  E, )  -  e  (sin  E2  -  sin  E, ) 

From  Equations  (F-4')  and  (F-10), 

rir2vce  sin2  °R  _  vce  _  "V 


cd 


F  ju2b2  juAF 

From  Equations  (F-4'),  (F-5),  and  (F-16), 


{sin  (E2  -  Ex)-e  (sin  E2  -  sin  E, )}.  (F-22) 


.  1  dF  a 2 

cotfcf, - - —  = - 

E  dOR  /- r2  sin  dh 


/to6_{— cos  ^  +  e  (cos  £’|  +  cos  E2 )  -  e2  (1  +  sin  E,  sin  E2 )} .  (F-23) 

ivE 


-{-  cos  (E2  -  E, )  +  e  (cos  E,  +  cos  E2 )  -  e2  (l  +  sin  E,  sin  E2 )} 
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From  Equations  (F-7),  (F-12),  and  (F-23) 
f  1  dF  A 


v 


COt  Or,  - 

v 


—  a  ^ {^2  “  ^ (sin F2- sinF,)}. 


r,  \  F 


3 

<2  V,. 


j-cos(F2  -  F,)  +  c(cosF,  +cos F2)  -c2  (1  +  sinF,  sinF,)} 

-F,  -c  (sin  F,  -sinF,)}{-cos(F,  -F,)  +  c(cosF,  +  cos  F,)-c2  (1  +  sinF,  sinF,)}. 
v  F 

(F-24) 

Substituting  Equations  (F-17),  (F-19),  (F-21),  (F-22),  and  (F-24)  into  Equation  (F-15), 

IF,  -F,  -c(sinF,  -sinF,)}j-cos(F,  -F,)  +  c(cosF,  +  cosF,)  -c2  (1  +  sinF,  sinF,)} 
+{sin(F2  -F,)-c(sinF,  -sinF,)}{l-cos(F2  -F,)-c(cosF,  +  cosF,)  +  2c2} 

+  sin(F2  -F,)|cos(F2  -F,)-c(cosF,  +  cos  F,)  +  c2} 

+  (F2  -F,){cos(F2  -F,)-c(cosF,  +  cosF,)  +  c2  (l-2sinF, sinF,)} 
sin(F2  -F,){ 3c2  +  l-2c(cosF,  +  cosF,)} 

-c(sinF,  -sinF,){l  -  2cos(F2  -  F,)  +  c2  (l  -sinF,  sinF,)}-3(F,  -  F,)c2  sinF,  sinF2 


\^dRjv 


_  QjTg 


Vf 


a3v 


Vf 


Now 


cosF,  +cosF2  =  2cosvcos(v  +  F,)  =  2cotxsinxcos(x  +  F,)  =  cot  v  (sinF,  -  sinF,) 
sin(F2  -F,)(cosF,  +cosF2)  =  2  sin  v  cos  v  cot  v(  sinF,  -sinF,)  =  2cos2  .v(sinF,  -sinF,) 

=  {l  +  cos(F2  -F,)}(sinF2  -  sinF,) 


S',  1 

_ 

-c(sinF,  -sinF,) (3  +  c2  -c2  sinF,  sinF, )  +  (3c2  +  l)sin(F,  —  F, ) 

Vf 

Vrfl 

-3c2  ( F,  -  F, )  sin  F,  sin  F, 

fdv^  dvJdtf^ 


c9 


\d9Rj, 


dv 


c8 


dt. 


\dORj 


which  yields  Equation  (F-14). 


dtf 


Equation  (F-14)  is  analytic  everywhere  provided  that  — —  ^0.  It  will  now  be  shown  that 

dv„ 


v  c8 


dt. 


dv 


<  0  universally.  From  page  18  of  the  main  text, 


C0 


dtf 

lim — —  = 

x^°  dvc0 


ry +  r2 -2  Jr\r2  cos - 


5jLlr2  sin  0K 


0, 


1 J  (ri  +  ri  )2  -  riF  cos2  y  +  (F  +  ri ) < cos y  \  <  0, 


0B 
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which  proves  the  proposition  for  a  parabolic  trajectory.  For  an  elliptical  trajectory,  since 

dtf  z\x  dv  a  dtf 

— —  =  ,  where  — —  <  0  (from  Equation  (E-l)),  it  is  desired  to  prove  that  — —  >  0.  Now, 

dvc0  dx  dx 

dx 

from  Equation  (5)  of  the  main  text, 


dtt 


<)X  yJljUS 


CSC  X 


iri  +  r2 )2  {3csc°  x (l  —  xcotx)  - 1} 

+2rjr2  cos2  — |6csc2  jc-5  -  3xcot  x(2csc2  x  - 1)| 


{rx  +  r2)cos^-|3xcscx(4csc2  x  -3)  -  (l2csc2  x  -ljcosxj 


To  prove  that  the  above  expression  is  positive,  the  following  theorem  will  be  employed:*  For  all 
real  values  of  u,  the  expression  Au  +  Bu  +  C  has  the  same  sign  as  A,  except  when  the  roots  of 

the  equation  Au  +  Bu  +  C  =  0  are  real  and  unequal,  and  u  has  a  value  lying  between  them. 
Putting 

u  =  1 . — l  sec  — ,  A  =  3csc2x(l-xcotx)-l,  B  =  3xcscx(4csc2 x -3) - (l2csc2 x - l)cosx, 


r\ri  2 


C  =  2|6csc2  x-5-3xcotx(2csc2  x-l)}. 


we  have 


dt 


f  _  rlr2 


cos2  —  esc  x[a.u2  +  Bu  +  c). 


dx  V2^  2 


Clearly,  u  >  1.  Since  A  >  0  (from  Equation  (19)  of  the  main  text),  the  above  theorem  implies 

dt  f  ~  2 

>0  when  either  iT-4AC<0,  or  5“-4AC>0  and  u>w,  where 


dx 


w 


1  m - 'l 

—  B  +  V B~  -  4 AC  J.  Performing  the  computations, 

2  A 


B 2  -4AC  =  9x2csc2x  +  6xcotx  +  cos2x-16  =  (3xcscx  +  cosx)  -16 
=  (3xcsc  x  +  cos  x  +  4)  (3xcsc  x  +  cos  x  -  4) . 


Now 


—  (3x  esc  x  +  cos  x  -  4)  =  3csc  x  (l  -  xcot  x)  -  sin  x  =  Asin  x  >  0,  0  <  x  <  n . 
dx 


*  Hall,  H.  S.  and  Knight,  S.  R.,  Higher  Algebra,  Fourth  Edition,  Macmillan,  London,  1936. 
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Thus,  3x  esc  x  +  cos  x-4  increases  monotonically  in  the  interval  [0, 7l) .  Since 
lim(3vcscv  +  cosv-4)  =0,  it  follows  that  3* esc  v  +  cos  v  -  4  >  0,  where  0  <  v  <  n .  Hence, 

x->0  v  7 

B 2  -  4AC  >  0.  By  expanding  cos  v,  esc  x,  and  cot  x  in  powers  of  x,  it  can  be  shown  that 

V5-1 


limw 

x— >0  2 


Also, 


f  12  —  sin2  v) cos  v sin  *-3.*  (4 -3 sin2  v)  +  sin2  xJ(3x  +  cos  vsin  x)“  -16  sin2  x 

lim  w  =  lim- - - - ^ - - - - =  -2  ■ 

2[3(sinv-vcosv)-sin  v} 


Plotting  w  versus  x  for  0  <  x<n  reveals  that  w  decreases  monotonically  in  the  interval  [0,  k)  . 
Hence, 

V5  —  1  dtf 

-2  <w< - <  1  <  u  =>  — —  >  0. 

2  dx 


dt 


f 


For  a  hyperbolic  trajectory, 


dtf  _  3y  „  _  I ^ 


- - .  Differentiating  vcg-- 

ovcg  ovC0  V  r{S  2 


f  sin  — ,  where 


dy 


s'  =  #j  +  r2  -  2^/^V cosh  y  cos  ,  we  have  ^  ,2 

dtf 

to  prove  that  — —  <  0.  From  page  16  of  the  main  text, 
dy 


r2  sinh  y  sin  0R  >  0.  Thus,  it  is  desired 


— —  =  ,  i  csch  y 
dy  J2/IS' 


(r,+r2)  {3csch2y(ycothy-l)-l} 

+2^r2  cos2  —  {3y  coth  y  (2csch2y  + 1)  -  6csch2y  -  5} 


+sfr\K,  (r]  +  r2)cos— |(l2csch2y  +  l)cosh  y  -  3y  csch  y  (4csch2y  +  3)j 


Vl 


e„ 


^  —cos2-^-cschy(A'u2  +B'u  +  C '), 


where 


A'  =  3csch2y  ( y  coth  y  - 1)  - 1,  B'  -  (l2csch2y  + 1) cosh  y  -  3_y  csch  y  (4csch2y  +  3) , 
C'  =  2|3y  coth  y  (2  csch2  y  + 1)  -  6  csch2  y  -5|. 


F-ll 


NSWCDD/TR-03/8 1 


Since  A' <  0  (from  Equation  (22)  of  the  main  text),  the  theorem  on  page  F-10  implies  that 

— —  <  0  when  either  B'2-4A'C'<0,  or  B'2-4A'C'>0  and  u>w\  where 

dy 

w'  =  -  4 B'2  -  4A'C,]j .  Performing  the  computations, 

B'1  -  4A'C'  =  9  y2csch2  y  +  6 y  coth  y  +  cosh2  y  — 16  =  (3y  csch  y  +  cosh  y)~  - 16 
=  (3 y  csch  y  +  cosh  y  +  4)  (3 y  csch  y  +  cosh  y  -  4) . 


Now 


d_ 

dy 


(3y  csch  y  +  cosh  y  —  4)  =  sinh  y  -  3csch  y  ( y  coth  y  —  1)  =  -A'sinh  y  >  0,  y  >  0. 


Thus,  3  j  csch  y  +cosh  y  -4  increases  monotonically  in  the  interval  [0,°°).  Since 
lim(3ycsch  y  +  cosh  y -4)  =  0,  it  follows  that  3y  csch  y  +  cosh  y -4  >  0,  where  y  >  0.  Hence, 

T— >0  V 

B'~  -  4 A'C'  >  0.  By  expanding  cosh  y,  csch  y,  and  coth  y  in  powers  of  y,  it  can  be  shown  that 


..  ,  45-1 

lim  w  = - 

y^o  2 


However, 


(l2  +  sinh2  yjcosh  y  sinh  y  -3y  (4  +  3 sinh2  y)  +  sinh2  y  J(3y  +  cosh  y  sinh  y)  - 16 sinh2  y 

lim  w  =  lim  - — - - - - - ^ — - - — - — — - : - : - — 

2 [sinh  y  -3(y  cosh  y-  sinh  y)| 


and  a  snag  is  encountered.  This  snag  can  be  circumvented  by  observing  that 

Q  _  Q 

S'  =  Tj  +  r2  -  2444  cosh  y  cos  =  444  cos  2  coshy)  >0  =>u>  2coshy, 

/ 

w  , 

which  suggests  finding  lim - .  Using  the  above  expression  for  w,  it  can  be  shown  that 

y->°°  cosh  y 

w'  w'  45  —  1  / 

lim - =  1.  Since  limcoshy  =  l,  it  follows  that  lim - = - .  Plotting  w  sech  y 

cosh  y  cosh  y  2 

versus  y  for  0  <  y  <  °°  reveals  that  v/sech  y  increases  monotonically  in  the  interval 

Hence, 

45—1  ,  ,  u  ^rf 

- <  w  sech  y  <  1  =>  w  <  cosh  y  <  —  <  u  — —  <  0. 

2  ’  2  dy 
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Thus,  it  has  been  proven  that  — —  <  0  universally.  A  heuristic  argument  that  this 

dvc0 

inequality  is  indeed  valid  was  presented  to  the  author  by  William  Davis  in  a  private 
communication.  If  either  vcr  or  vc9  is  increased  while  holding  the  other  constant,  dR  will 

increase.  Therefore,  to  keep  0R  constant,  either  vcr  has  to  be  increased  while  decreasing  vc0  or 
vice  versa.  However,  simultaneously  increasing  vcr  and  decreasing  vc9  tends  to  increase  !  f , 
while  simultaneously  decreasing  vcr  and  increasing  vcQ  tends  to  decrease  l  f .  In  either  case, 
dtf 


To  recapitulate,  given  rl5  r2,  and  vcQ  ,  the  desired  algorithm  is 

f  r  \  u  0  A  r  A  Q 

eR  =  arccos (r,  •  r, ) ,  v„  =  cot  0R  — L esc 0R  vc8  + - tan ,  vtr  =  1  — L  vc0  cot -*-■ -  vcr , 

{  r2  J  Wee  2  l  ri)  2 

^  =  j2^snr|- rx  v;e  | 'n  +  r2-  2 ^  cos  1 1  rxv]0  \  rx  +  r2  +  2^_  cos  ^  -  2  jur2  sin2  ^  j , 

_  Hr\rlvlt>  s^2  Or  h_Kr2v2c8sm9R 

CL  9  Cf  I -  ,  C  1  9  9 


esinE.  =  r^2L,  esin£^  =  -5^=,  E~  —  E.  =  2arctan- 


( ri  +  r2 )  vle  ~  2Ar2  sin2 


(^t  fl 

— —  =  °-j=  {sin  [E2  -  Ex)  +  e^ini?,,  -sin£'1)}{3(£’2  -  Ex)-e(sinE2  -  sin  .E,)}- 16  sin2 
dvc0  VT  L 

a 3  3v  „  \e  (sin  E2  -  sin  iq )  (3  +  e2  -  e2  sin  Ex  sin  E2 ) 

— ~  =  —j=  vc-e  — ~  )  f  • 

v d0R  Jif  vE  dtf  -(3e2  +  l)sin(E2 -E{)  +  3e2(E2  -  E[ )  sin  E[  sin  E2 


^2  Ei 
2 
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APPENDIX  G 

A  GENERAL  EXPRESSION  FOR  THE  NULL-MISS  DIRECTION 
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The  null-miss  direction  is  defined  by 


where 


P/w  - 


d\c 

dt( 


\d\c 

dtf 


d\  dv 


dtf  dtj 

dv,  dv 


-  + 


dv 


dr. 


(K 

dt . 


-r,  + 

dtf  dt 


dv 


c0 


V 


dt 


0p 


7 


_  _^r  =  Q  x  r, , 

dtt  2 


^dv  ^ 


vdr2y 


vdr2  ; 


7 

^dv_  A 


v^y 


A 


0 


dv 


C0 


v  ar2  y 


l 


dv 


y 


c6> 


dtf„ 


'2  vucw 
o 


r2  sin  0R 


cos  9r  sin  9r  0 
-  sin  9 R  cos  9 R  0 
0  0  1 


(G-l) 


ce A 

\^9r  J 

V  dr2  J 


is  given  by  either  Equation  (E-16)  or  (F-14), 


fog 

k^9r  y 


is  given  by  Equation  (E-17), 


is  given  by  Equation  (E-19),  and 


<K 

V  dr2  J 


is  given  by  Equation  (E-20).  We  will  now 


d\ 

show  that  — -  =  -Avc ,  where  Avc  is  given  by  either  Equation  (7)  or  (8)  of  the  main  text. 


dt 


Suppose  that  the  target  vector  r2  is  perturbed  by  a  small  amount  Sf2  (Figure  G-l).  Let  r2 
be  the  new  target  vector,  89 R  =  9'R  -  9 R ,  and  SjB  =  jB'  -  jB .  Then  Sr2  =  r2  -  r2  ~  r2S9RQ2  +  8rp n , 
where  sgn  [Srp )  =  - sgn  {SfB) ,  the  minus  sign  being  attributed  to  the  fact  that  a  positive  change  i 
! 3  points  in  the  opposite  direction  as  n .  From  Figure  G-l,  we  have 


in 
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Figure  G-1 .  Representation  of  a  small  perturbation 
in  the  target  vector 


Sr 

sin — —  =  -sin#'  sin(£/?) . 

r2 

If  these  finite  changes  are  made  infinitesimal,  then 
Sr  dr  dr 

sin — —  =  sin — —  =  — — ,  sin($?)  =  sin d j3  =  dj3,  0'R  =  dR  +  dOR  =  0R 


dr 

dr  =  -r2  sin  0Rd  jd.  or  — —  =  -r2  sin  0R. 

djB~ 


From  Equation  (G-1),  we  have 


dvc  _  v 


cd 


drp  r2  sin  0R 


n 


d\c  _  dvc  drp 


djB  dr  d  (5 


~Vc9n  ■ 


We  will  now  evaluate  — ^ -  and  — —  for  which  we  need  the  following  identities:’* 

dtf  dtf 

cos  0R  =  sin  cp  sin  cpT  +  cos  cp  cos  cpT  cos  AD , 
sin  0R  sin  j3  =  cos  (pT  sin  AD , 
cos  (pT  cos  Ad  =cos(pcos0R-  sin  cp  sin  dR  cos  j3, 
cos  cpT  cos  j3b  =  sin  <p  sin  dR  -  cos  ^cos  0R  cos  j3, 
sin  cpT  =  sin  (p  cos  0R  +  cos  <p  sin  0R  cos  /?, 


(G-2) 


(G-3) 

(G-4) 

(G-5) 

(G-6) 

(G-7) 


Sofair,  Isaac,  K40  Training  Guide  4085.1,  A  Detailed  Derivation  of  Formulae  Arising  in  Spherical  Trigonometry, 
Naval  Surface  Weapons  Center,  Dahlgren,  VA,  May  1986  (currently  Naval  Surface  Warfare  Center,  Dahlgren 
Division). 
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where  AD=AT-A,  X  and  Xr  being  the  longitudes  at  release  and  impact,  respectively;  (p  and  <pT 

are  the  corresponding  spherical  latitudes  (i.e.,  latitudes  corresponding  to  a  spherical  earth),  and 
fdh  is  the  back  bearing.  Differentiating  Equation  (G-3)  with  respect  to  tf  and  using 

d/L 

Equation  (G-4)  and  the  fact  that  — —  =  Cl .  we  have 

dtf 


■  a  dOR 
sin 


dt. 


■  -  cos  <p  cos  cpT  sin  XD 
dd , 


dX 


D 


dt. 


-Cl  cos  cp  sin  0R  sin  fd 


dtf 


Cl  cos  (p  sin  fd . 


(G-8) 


Differentiating  Equation  (G-4)  with  respect  to  tf  and  using  Equations  (G-5),  (G-6),  and  (G-8), 
we  have 


=>  sin  0R  cos  (d 


d/d 

dtf 


=  Cl  (cos  cpT  cos  XD  -  cos  ^cos  0R  sin2  / 3 ) 

=  Cl  (cos  <pT  cos  XD  -  cos  (3  cos  0R  +  cos  ^cos  0 R  cos2  y 3 ) 
=  Cl  (cos  <p cos  0R  cos2  fd  -  sin  <p sin  0 R  cos  fd ) 

=  -Cl  cos  (pTcosfd  cos  /dh 

esc  0R  cos  (pT  cos  fdh . 

dtf 


(G-9) 


Now  if  P  is  a  unit  vector  in  the  direction  of  the  North  Pole,  then  =  PG  . 

P  x  r2  •  n  =  P  •  r2  x  n  =  -P  •  02  =  P  •  (rt  esc  dR  -  r2  cot  0R )  =  P  •  rt  esc  dR  -  P  ■  r2  cot  0R 
=  sin  epese  0 R  -  sin  cpT  cot  0R 

=  sin  cp  esc  0R  -  ( sin  cp  cos  0R  +  cos  (p  sin  0R  cos  fd)cotOR  (from  Equation  (G-7)) 
=  sin  cp  sin  0R  -  cos  cp  cos  0 R  cos  fd  =  cos  <pT  cos  fdb  (from  Equation  (G-6)) 


d  fd 


Cl  x  i%  ■  n  =  Gr2P  x  r2  ■  n  =  —r2  sin  0R - (from  Equation  (G-9)) . 

dtt 


(G-10) 


Since  fd  is  the  angle  between  P  x  r,  and  r,  x  r2 ,  it  follows  that 


(Pxf1)x(f1xf2) 


Pxr, 


|r,  x  r2 1  sin  fd  =  cos  tp  sin  0 R  sin  fd 
|(Pxf,.f!)f,-(Pxf1.f,)f2|  =  |(Pr,xr2)r,-(p.r1xrl)f2 
P  •  x  r,  =  P  •  n  sin  dR  =>  P  ■  n  =  cos  cp  sin  fd. 
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Referring  to  page  7  of  the  main  text  and  Equation  (G-8),  we  have 


A0R  =  il  ■  n  =  -£2P  ■  n  =  -flcos^sin  j3  =  - 


dOR 
dl , 


(G-ll) 


(  dt  ^ 


1  + 


Ax  =  —- 


K^Or  j 


A  0R 


dx 


fdt  ^ 


dt 


f 

dx 


dt 


1- 


/ 


\d9Rj 


dOR 
dt  f 


(G-12) 


From  page  8  of  the  main  text  and  Equations  (G-10)  and  (G-ll), 

dr, 


dt  r 


■  Q x r2  =  r2  sin ^A^r,  -  r2  cos dRAdJbx  +  x r2  ■  n) n 


.  „  d  6>  n  dOR  -  .  _  d/3  „ 

■  -r,  smOD  — — ! -r,  +  r,  cos8D—JLQl  -  r2  sin#,,  — —  n  . 


'R  i.  *1  1  '2' 

dt s 


dt. 


dt , 


(G-13) 


From  Equations  (G-l),  (G-2),  and  (G-13),  it  is  easy  to  see  that 


^dv  ^ 


dr. 


A 


V  Jt/  dtf 


dv 


\dOR  J  dtf 


d  6„ 


-r,+ 


rx,.  \ 


cO 


dv 


V  Jt,  “7 


d6R  n  .. 

- o,  —  - n 

dtf  dt, 


(  X-.  \ 


dv,  dv„ 


dtf  dt 


-  + 


dv 


dr,  dv. 


/ 

f  x^-,  \ 


dv 


dOR  dvc  d  fd 


■  + 


yddR  Jt  dtf  djB  dtf 


f  V  ^*2  J,f  “7 


dt,  dt 


-  + 


dv„ 


dOR  dve  d/3 


f  “7 


■  + 


dt,  d/3  dt , 


dvcrX  I  v 

- r,  h - o,  —  vrfl - n  . 

dtf  dtf  dtf 


(G-14) 


Now 

dv„.  dv. 


dtf  dtf 


-  + 


f  dv.  A 


dO, 


d  0„  dv„.  dx 


R  _  XZcr 


\  wlyR  J,  ul  f 


dtf 


dx  dt. 


■  + 


rdv„  A 


\dQRj 


+  - 


dv. 


f  dx  x 


dx 


\ddRj, 


dOR 

dtf 


dvcr 

dx 

dv, 


dx 

dtf 


■  + 


f  dx  x 


\dOR)  dtf 


d  0„ 


+ 


^dv„,  ^ 


j 


dOR  _  dvcr  dx 
dt. 


f  dt,  ^ 


dx  dtf 


\dOR  j 


d  Or 
dtf 


•  + 


^  dv,  A 


\dOR  J 


dO R 

dtf 


fx,.  \ 


dx 


-Ax  - 


dv 


\ddRJx 


A 0R  (from  Equations  (G-ll)  and  (G-12)). 


Substituting  for 


dv 

cr 

dx 


and 


(dv^ 


j 


from  Equations  (E-4)  and  (E-5)  and  using  Equation  (8)  of  the 


main  text,  we  obtain 
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<K 

dt. 


2/i  . 

Vsinx 
S3 


6, 


i]  +  r2  sin  -  ^Jrtr2  cos  x cos 


+ 


\JH Lsin^ 

'2rxS3  2 


/,  sin  x  +  r2  -  Jrj r2  cos  x  cos  - 


Ay 


A#„ 


=  -Avc  ■  rj . 


(G-15) 


Similarly,  from  Equations  (E-l)  and  (E-3),  we  have 


**.  — *?*.  Ax- 


( :v.  A 


C0 


dt. 


dx 


dv 


k^0r  j 


AO, 


JL 

2S3 


n  sin  x  sin  dDAx  ■ 


Mri 
2  rxS3 


t  \  Or  I - 

(rt  +  r2  j  cos  -  ^/r,r2  cos  x 


l  +  cos°  — 


AO,., 


J 


-Avc  ■  0, . 


(G-16) 


From  Equation  (F-7)  and  the  expression  for  g  in  Appendix  C,  namely, 


i]  r2  . 

g  =  ~ p=  sm  0R  , 


4i*p 


we  have 


8  = 


_  r2  sin  0R 


(G-17) 


vc9 


Substituting  Equation  (G-17)  into  Equation  (G-10),  we  obtain 


Q  x  r,  n  dB  „ 

- 1 —  =  -vrf  —  =  -Av„  •  n 


8 


dt , 


dp 


^vce  —  =  Av„-h. 


dt 


(G-18) 


/ 


Substituting  Equations  (G-15),  (G-16),  and  (G-18)  into  Equation  (G-14),  we  finally  obtain 

d\. 


dt 


f 


-(Ava,r1)f1-(AvI.0,)e,-(Ayc-n)n  =  -AvI, 


which  proves  the  identity. 

Note  that  the  expression  for  on  page  G-3  is  general  in  that  it  is  applicable  to  all  three 
types  of  trajectories. 
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APPENDIX  H 

DETERMINATION  OF  FOOTPRINT  ACHIEVABILITY  IN  CONNECTION  WITH 
TIME  OF  FLIGHT  AND  CANTING  ANGLE 
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Assume  a  scenario  consisting  of  two  targets  whose  spacing  is  8s .  Let  r,  be  the  position 
vector  at  first  release,  r2  the  first  target  vector,  r2  the  second  target  vector,  0R  and  [3  the  range 
and  bearing  corresponding  to  the  first  target,  and  0'R  and  f3'  the  range  and  bearing  corresponding 
to  the  second  target  (Figure  H-l). 


Figure  H-1 .  Canting  angle  as  a  measure  of  the 

deviation  of  r'  from  the  downrange  arc,  0 


The  canting  angle,  ac ,  measured  clockwise  from  the  downrange  direction  at  first  impact,  is 
the  angle  between  the  trajectory  plane  formed  by  r,  and  r2  and  the  plane  formed  by  r2  and  r2 . 
Its  domain  is  0  <  ac  <  In .  It  is  easy  to  show  that 

sgn($?)|(?,  x?2)x(?2  xr2)|  sin 0'R  sin 

a  =arctan - - — 1 — — - - - L  =  arctan - - - - — ,  ,  . 

(r,xf2)-(f2xr2)  cos  0 R  sin  0 R  cos  ( 8J3 )  sin  0 R  cos  0 R 

For  a  given  canting  angle,  ac ,  it  is  possible  to  predict  the  relative  capability  of  a  missile  for 
different  times  of  flight  by  determining  the  relative  sizes  of  the  velocity  gains  necessary  to 
produce  a  given  increment  in  impact  position.  Minimum  velocity  gain  reflects  optimum  missile 
capability.  Specifically,  the  magnitudes  of  the  minimum  velocity  changes  necessary  for  the 
given  times  of  flight  must  be  compared.  The  minimum  velocity  change  is  achieved  by 
considering  the  effect  of  changing  tf . 

Suppose  that  a  definite  trajectory  defined  by  rp  r2,  and  tf  has  been  selected,  as  well  as  a 
fixed  small  change  in  impact  position,  Sr2 .  For  any  small  change  8tf  in  tf ,  we  have,  using  the 
results  of  Appendix  G, 
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d\  a 

£>v  =  —^8t  +  —^Sf  =-P 

U'c  U‘f  -i-  Ul2  *DK 

dtf  dr, 


dv 


dtf 


<K 

dp 


dv  dv 

where  it  is  implied  that  — -  and  — -  are  evaluated  at  constant  tf .  The  value  of  St  f 
V  dOR  dp  f  f 

_  _  dv 

minimize  Svc  is  the  one  for  which  S v  is  perpendicular  to  — i.e., 

dt. 


Svc-^  =  0,ov  S\c-PDK=0  => St 

dtf 


dyc 

-v 

dtf 

V 

Vdk-^Mr+Vdk-^SP 

DK  doR  R  DK  dp  ) 


Substituting  this  value  of  Stf  into  the  expression  for  Sv c ,  we  obtain 


(<?vc)  .  = 

v  L  '  min 


<K 

d  0„ 


p 

1  DK 


dv,  ^ 


dO, 


SO,,  + 


r  J 


fog 

dp 


*  A 


p 

1  DK 


dv 

Jp 


■sp. 


If  these  finite  changes  are  made  infinitesimal,  then 

d(vc)  dv  A  f 

dOB  d0„  DK 


dv 


A 


d(vc)  dv  A 

v  c  '  min  c_ p 


( 


dp  dp 


DK 


d\ 


Itt  A 


DK 


dp 


=  vce  {pd*  (PDr  ■  n)  -  n},  using  Equation  (G-2). 


From  Figure  H-l,  we  have 
Sr, 


sin — —  =  -sin#'  sin(£/?)  =  -sin  —  sinac,  cos — -sm(SOR)  =  sin— cosac . 


ri 


r2 


r2 


r2 


If  these  finite  changes  are  made  infinitesimal,  the  above  equations  reduce  to 

_  Mr  _  cos  ac  dp  _  sin  ac 


dr  =~r2  sin  0Rdp  =  -ds  sin  ac ,  r2d 0R  =  ds  cos  a 


ds 


r2  ds  r2  sin  0, 


Hence, 


d(vc)  .  d(vc)  .  dOR  d(vc)  dp  l\ 

v  c  '  min  v  c  '  min _ R_  _j _ v  c  '  min  /  J 

ds  dOR  ds  dp  ds  r2  I 


d(vc)  .  sin  a  d(vc) 

v  c '  min  _|_  c  v  c '  r 


1 


cosa. 


5v c  p 

dOB  DK 


0v 


1  DK 


v 


sin  a, 

■  H - —  V 

sin  0P 


dOR  sin  0R  dP 

{P«(P«n)-n} 


H-4 


that  will 
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With  the  knowledge  that  a  lower  value  of  -  111111  implies  a  greater  missile  capability, 

ds 

d  (v  )  . 

the  values  of  - —  corresponding  to  a  specific  canting  angle  for  different  times  of  flight  can 

ds 

d(v  )  . 

be  compared.  Thus,  - —  can  be  used  to  determine  the  dependence  of  missile  capability  on 

ds 

the  time  of  flight. 
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APPENDIX  I 

COMPUTATION  OF  THE  RATE  OF  CHANGE  OF  THE  NULL-MISS  VECTOR 


I- 1/1-2 
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dPn 


ell \ 


Since  PDK  is  the  direction  in  which  the  guidance  system  is  attempting  to  steer  the  missile, 
-  predicts  how  the  missile  is  turning,  i.e.,  how  the  direction  of  thrust  is  changing.  In  this 


appendix,  a  closed-form  expression  for 


dP„ 


dt, 


in  the  local  coordinate  frame  will  be  derived.  The 


results  of  Appendices  E,  G,  J,  and  K  will  be  used.  In  Appendix  G,  it  was  shown  that 
-  d\„ 


L  DK 


dtt 


f^cr  +dVcr  cWR ' 


V  dtf 


ddR  dtf  j 


ri~ 


fdvce  +dvc6  clOR 


Kdtf 


dOR  dtf  j 


0i+^  — n 

dtf 


(1-1) 


where  anc]  f}h±  impiy 


d  0„ 


d6v 


' R  w  R 

respect  to  tf  can  be  expressed  as 


fdv  ) 

f  dv.^ 

and 

d0Rj 

tf 

Kd0R; 

respectively.  The  derivative  of  Pj)K  with 


dP, 


DK 


dtt 


dVDK  +  d PDK  elr[  _  d P/M.  +  dP/)A.  ih\  +  dPDK  dOR  +  dP^  el/3b 


dtf 

3Pn 


dtf 


-  +  v„ 


dr,  dtf 

3P, 


dt. 


dr,  dt. 


dOR  dtf 


d Ph  dtf 


-  + 


dPr>K  ddR  ,  dPnK  d[l 


+  - 


dr,  dOR  dtf  d/3h  dt. 


(1-2) 


where  the  partials  are  with  respect  to  a  fixed  target  r2  and  a  changing  initial  position  r, . 


dP  dP  dP  dP 

We  now  proceed  to  derive  expressions  for  DK  ,  DK  ,  — — ,  and  — —  .  The  first  two 

dtf.  dr,  dOR  dj3h 

partials  are  straightforward.  A  change  in  tf  without  changing  r,  does  not  change  the  local 
coordinate  frame.  Since  f,  and  r2  are  both  constant  in  this  case,  0 R  and / 3  -  and  therefore 
clOR  clB 

and  — - remain  constant 


d  dOR  _  d  d/5 


dtt 


dtf 


dtf  dtf 


dtf  dtt 


0.  Hence, 


dP 


DK  _ 


dt. 


V  dtf 


-  +  - 


d\r  d0R 

dORdtf  dtf  j 


cQ 


V  dt) 


!  d2vce  dOR  (  dvcg  d/3  ^ 


dORdtf  dtf  j 


dtf  dtf 


(1-3) 


Similarly,  a  change  in  t\  does  not  change  the  local  frame.  Since  0R  and  fd  -  and  hence  their 


derivatives  -  are  independent  of  r, ,  ^  C^R 

dr,  dtf 


d  clfd 
d  r,  dtf 


■  0.  Hence, 


dP 


DK  _ 


dr, 


dS 


■  +  ■ 


d2v„„  ddD 


drxdt f  drxdOR  dt 


f  J 


c6 


d2v  g  clOR  . 

+  r—rsi - -  10,  + 


drxdtf  drxdOR  dt 


f  J 


dvcg  dP 

dr,  dtf 


n . 


(1-4) 
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3P 

In  computing  — — ,  observe  that  a  change  80  K  in  0R  is  due  to  a  change  in  r,  so  that  the 

deR 

local  coordinate  frame  changes.  Let  r',  0',  n  be  the  new  local  coordinate  frame.  Then, 
r'  ~  r,  cos  ( 80  R )  -  0,  sin  (80 R ) ,  the  minus  sign  being  attributed  to  the  fact  that  a  positive  change  in 
0R  points  in  the  opposite  direction  as  0, .  If  80 R  is  made  infinitesimal,  then 
f '  =  r j  —  0 j cl0R,  0'  =  nxr' =  0,  +rxd0R  .  Let  P'DK  be  the  counterpart  of  P/)A,  in  the  new  frame. 
Neglecting  terms  in  ( d0R)~ ,  we  have 


3P 

p'  p  i  17  DK  .in 

rDK  rDK  '  -  UUD 


d0D 


\(Pok)  + 


H Pqk l 
d6„ 


f  P  +  I  (  PDK  L  + 


d  (PDK)6 
d0„ 


dOR  >0'  +  \{PDK)n  + 


ai^detyn 
d0,, 


(/’«),  +  ]  (f,  -  w) + |(p„ ) 


.  d  ( PDK  )6 
e  deD 


dOR  H  0,  +  rxdOR} 


+\(P»K),  +  dP^JLde.',ii 


?DK  + 


d(PDK)r 

d0„ 


+  {PDK)g  kl  + 


d  (P/  jk  )6 
d0„ 


-(PDK)r  9i  + 


d(PDK)r 

d0„ 


n 


cWD 


ap 


DK  _ 


a  0„ 


d0„ 


+  ( Pdk)q  f  T  + 


Ja  (/>«■), 

<  - - 

d0R 


(PDK)r\^  + 


d  (Pqk  )j 


n 


d\ 


-  +  ■ 


dv 


(  ^2. 


c8 


dORdtf  dtf 


+ 


■  +  ■ 


dv 


c0 


V  d0R  66 R  J 


d0D 


dt, 


d2vce  dvc 


deRdtf  dtf 


■+ 


rd2 


vce  dvc 


\  d0]  00 R  J 


d  0„ 


dt 


0i 


/ 


+ 


rdvstLd_P_  +  v  a  d_P' 
d0D  dtf  3  0,,  dt 


^  yj  ys  ft  mi- y 


n. 


'R  w,7  J 


(1-5) 


Note  that  Equation  (1-5)  should  contain  additional  terms  involving  ;  however,  it  is 


dOR  dtf 


shown  in  Appendix  K  that  =  q 

dOR  dtf 


1-4 
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Similarly,  in  computing  — — ,  a  change  in  [3b  is  due  to  a  change  in  r, .  Let  r',  0',  n'  be 

°Pb 

the  new  local  coordinate  frame.  Consider  a  perturbation  £r,  resulting  in  a  small  change  8(3 b 
without  changing  0R .  By  an  analysis  similar  to  that  on  pages  G-3  and  G-4  for  a  perturbation  in 
the  target  vector,  d'r,  =  r'-Fj  ~  Srlpn  .  Since  a  positive  change  in  (3b  points  in  the  same  direction 

Sv  /v  /v 

as  n,  sgn(£r1/;)  =  sgn(£/?fc)  sin — —  =  sin#R  sin(^/?i) .  Also,  0'2  ~  Q2cos(S(3h)  -  nsin(£/?fo) . 

ri 

If  these  finite  changes  are  made  infinitesimal,  then 

drlp  =  /,  sin  dRd (3h  =>  r '  =  r,  +  n  sin  0Rd  (3h , 

0',  =  02  -  n d (3h  =>  0'  cos  0R  -  r'sin  0 R  —  0,  cos  0R  -  r,  sin  0R  -  n d (3b 

=>  0'  cos  0R  =  0,  cos  0R  +  (r' -  fj )  sin  0R  -  n d (3b  =  0j  cos  0R  -  n  (l  -  sin2  0R )  d (3b 
=  0,  cos  0R  -  n  cos2  0Rd (3h  =>  0,  =  0,  -  n  cos  0Rd (3b , 

n^r'xOj  =(fj  +  nsin#R<i/?fo)x(01  -  ncos0Rd  (3b^  =  n  +  (0,  cos^  -  r,  sin  0R )  d  (3h  =n  +  d2d/3h. 

Again  letting  Y'DK  be  the  counterpart  of  YDK  in  the  new  frame  and  neglecting  terms  in  (d J3b)~ , 
we  have 


-  -  dP 

p'  =p  +  DK  d  B 

A  DK  lDK  ^  -s  n  u  Pb 

°Ph 


3A 


9  (  PDK  )f, 

Wb 


teJ k+  (Pdk)  +-^i±dfib  e;  +  (pDK)+—^d/3b  U' 


d(PDK)f 

Wb 


d(PDK)r 


d(  PDK  )f/ 


( Pdk)t+  r  d (3 h  (fj  +  n sin 0Rd /?„)+( PDK )e  +  DK  8  d/3M-n cos 0Rd (3„ ) 


+  1  {Pdk)„  + 


d{PDK)r 


d(3b  i|n  +  (0jCos^  -fjSin 0R^d(3^ 


-Pdk  + 


d(PDK), 


-{PDk)s™9r  ki  + 


d(PDK)b 


+  {Pdk)„c™0r  01 


d(PDK)r 


+  {Pdk )r  sin  6r  ~{PDK )  cos  0R  n 
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dP 


DK  _ 


dA 


^  (  Aw  )r  (  n  \  •  1  *  _L  ^  (  Aw)#  ,  (  n  \  /J  1  n 

“  ( Aw  )  „  Sln  ®r  f  ri  1  *"  (  Pdk  )  COS  f  0i 


3A 


+  • 


3A 

^{Pdk)  j 

~W~ 


+  {pDK)r  sin 0R  ~(PDK )  cos eR  n. 


Since  vcr  and  vc6  -  and  hence  their  derivatives  -  do  not  change  with  either  [3  or  [3h  when  0R  is 
fixed,  it  follows  that 


dv 


c6 


d  dvcr  _  d  dvcr  _  d  dvc6  _  d  dvcff  _  Q 


dA  d/J*  c)[3h  c)0r  c)Ph  dtf  d/3h  d0R 


Hence, 


dP, 


DK 


dA 


+  iv 


dv£^J_d0JL+  in0  dfi 
yd0R  d/3h  dtf  ce  Rdtfy 

d  d  (3  dvc0  „  dv. 


ri- 


'«kl4_v  ms0  11 

v3ft3ftA,  '  dt. 


cG 


-x  ^  +  — — cos#„  ^sin#.  + 

d/3h  dt ,  dt ,  f)/  , 


dv. 


dv.. 


— ^cos#„  -^^sin^ 
vd6>s  d6>R 


d  0„ 


(1-6) 


dt. 


n. 


Substituting  Equations  (1-3),  (1-4),  (1-5),  and  (1-6)  into  Equation  (1-2),  we  finally  obtain 


dP, 


DK  _ 


dt. 


d2v  d2v  L  d2v 


dt2 


-  + v„ 


dr/dt. 


-+  2: 


r)0Rdtf 


-  +  v„ 


d2r 


-  +  - 


dv 


f  ^2. 


c0 


drjd^  d  t 


+ 


f 


-  +  - 


dv 


C0 


V  d  6>a,  , 


d0R  I  d0R 


dtf  dtf 


+ 


dv  d  d 0S  d/3 

— ^ - -  +  v„  sin  0, 

cB  R  dtf 


d  0R  d/3b  dt f 


3ft 


dt. 


d2! 


c6 


dt) 


+  V„ 


vc9 


dr/dt. 


+  2 


d2! 


c6 


d0Rdtf 


+  v„ 


d2E,  dv. 


drjd^  dt 


-  + 


^d2vc*  dv.. 


/ 


V  dd2  d0R  j 


d0R  d0R 


dtf  dtf 


+ 


3v^_3_rfft  3ft 

v  3ft  3ft  dt,  "  dt, 


3ft 


r/t. 


0, 


+ 


V6> 


d  d/3  +  d/3h  d  d/3 
ydtf  d  0R  dtf  dtf  d[3h  dt f 


+ 


dv  e  ,  dv,,e  d0R 


^  dtf  d0R  dtf 


+  v_ 


dv 


eft 


dr, 


d/3 


>  y 


dt. 


+  < 


dv 


dv. 


-^-cos0R  -- 
dt,  R  dt. 


■  si  n  0R  + 


dv 


dv. 


—  COS0R  -- — 

d0D  d0„ 


sin  0,, 


V  R 


d0R  d/3h 


"f  s 

where  all  the  required  derivatives  are  derived  in  Appendices  E,  G,  J,  and  K 


dtf  dtf 


n. 
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APPENDIX  J 

DERIVATION  OF  SECOND  PARTIALS  OF  vc*  AND  vcr 
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Put  vcr  =  £(vc0(tf,0R,rl'),0R,r1).  Using  the  equations 


f  r  > 

cot  6^  -  —  csc#s 

V<*  + 

A 

tan^, 

v„.  = 

r 

1- 

X 

a 

l  E  j 

2 

V 

r2y 

vc6cot- 


0 B 


we  have 


dvc0  vce 


1  +  ^ 

V  r2j 


0  dC 

tan  —  ,  ■z-^-  =  -v csc0R  (from  Equation  (E- 17)), 


2  ’  d0D 


dv„ 


dv 


c8 


f  r  A 

1-1 

V  hj 


d2C  _  1 

dKs  Vc0 


dC 

_  Vc,  +  Vtr  dC 

V 

—  cr 

f 

1  +  - 

V 

dvce 

Vcg  dvc0 

V  } 

tr 

d2c  _ 

dvtr  Q 

—  "  nen  h-i 

dvc8 

Vc6) 

dvcgd0R 

dvc8 

ri 


0R 

tan  —  , 


d2c 

d0l 


n  dC  1 
v.  cot  0R  +  — —  +  - 

R  d0D  2 


(  r ' 

1-1 

V  rU 


vcScsc2^-\csc0r 


=  iv 


1  6  0  ] 

.  (cot  0R  —  esc  0R )  +  -  ( v„  +  vlr )  tan ^csc2  yj  esc  0R 


=  R  (cot 0r-csc0r)  +  {vcr  +  vtr ) esc 0R }csc 0R 

=  (  Vcr  +  Vtr  COS  0R  )  CSC  2  0R  . 

dvcr  _  dC  dvc0  dv,r  _  dvrr  dvc0 


dtf  dvc0  dtf  ’  dtf  dvc0  dtf  ’ 


dv. 


d0R  d0R  dvc0  d0R  ’  dr, 


df.+  df  dvrt  dvc 


+  vc0  cot  0R 
dv..a  d0  cd  R 


1 


(Equations  (E-17)  and  (E-18)), 


dv. 


dv. 


f 


tj_  _  r  vc8  _  VL 

•  1 


d0D  1  dr. 


y2Vcd  +  d0R 


dv...  3  dv 


A 


2  1  dt 


c8 


(by  Equations  (E-13)  and  (E-19)). 


/  J 


Another  expression  for 


<K 

d0„ 


is  obtained  by  direct  differentiation  of  the  expression  for  vtr  above. 


We  have 


dv„ 


dv 


1 


d  0R  d  0 R  2 

d0R  2 


l-l 

V  r2j 

f  r\ 

1~1 

V  r2j 


2  0R 

V  0  CSC  — —  + 
cd  2 


1-1 

V  r2j 


dvn0  0,, 

— —cot  — — 

d0D  2 


2  @R 

v  0  esc’  — —  + 
00  2 


f  r  ^ 

1-1 

V  r2j 


cot  — 


0 R  dC  dv , 


cO 


2  dvc0  |  d0t 


=  -vcr  CSC  0R  +  ^V,r 

dvc0  d0R 
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d2tf 


dv 


is  obtained  by  differentiating  Equation  (14)  of  the  main  text.  To  facilitate  the  computation, 


cd 


it  is  convenient  to  define  some  recurring  terms.  Put 


Q  Q 

A  =  rl(rl+r2 ) ,  B  =  2rx  cos  -y-,  C  =  2  flr2  sin2  -y-,  D  =  4  jUtfr2  sin3  dR . 


2  2 

Then  F  =  {C  -  ( A  -  B)  v20}{(A  +  B)  v2ce  -  C}  and  Equations  (13)  and  (14)  of  the  main  text  can  be 
expressed  as 

Dvc0 


tf  = - 

f  IF 


|  AC-(A2-B2)v],  vf  Vf 

— - — - - - h— p^arctan 


dt 


7_  _ 


b2c 

-2  f 


f4f  Av-,-C 


H - i^aivuui - - - 

•If  av-c,-c 
■If  ' 


avce  r  ^vc 

Differentiating  the  second  expression  yields 

I"  «  [(A2-52)3v^-3C2(A2-52)(2A2-752)ve4e 

+8AC3  (A2  -  IB2 )  v2*  -  C4  (3A2  +  8  B2 ) 


2  CD  2 

r-2  Vc0 ' 


d2tf_ 

DvcS 

Ke 

Ui 

1 - 

B'C  [  +8 AC3  (A2  -252)v2^  - 

+3VF  |2  ( A2  -  B2 )  v4*  +  5 ACv)e  +  C 


,o  1052C2v4J  f  Vf 

- —  >  arctan  — - - 

F  J  Av^-C 


d2vc0  _  a  dvcg  _  a 

fV 

a 

r«'/i 

1  dvce  _ 

3d\ 

dtj  dtf  dtf  dtf 

dVc0 

l  ^Vc0  ) 

dtf 

1  dtf  ) 

dvc0 

d-vcd  _  d  dvc8  _  a 
dt2f  dt f  dtf  dt  f  ^  3vc 

Differentiating  Equation  (E-16),  namely, 

1 


c)v 

d0„ 


C0 


2  VcrVc0  ' 


r,  —  v 

1  cr 


r  r2  3  ^ 
— - tf 

Vv»  2  '  j 


dv 


cO 


d  t , 


cd 


\  riVcr ^ 

V  E ynj 


with  respect  to  t ,  and  0R  ,  respectively,  we  obtain 


a2! 


C0  _ 


dORdtf 


d\ 


dtf  1  dv 


C0  "V 


C0 


V  dt)  ^vcg  vc6  dtf  j 


dv 


c0 


dOD 


if 


+  - 


l  rAcr 


Wrr 


1  ~dVcr  rAcr  dVtr 


y  v cr  'dt,  r2v2r  dtf  j 


dv 


c0 


+  - 


d0R  Vc0 


ri  'Ch’cr  r2Vcr 


V  Vcr  &f 


v. 


2  dt 


-  v„ 


dv 


c0 


tr  ”‘.f  J 

a2 


dt 


i 


Vc0  dtf 


2Vcr  dtj  dv 


ce 
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d’l 


cO  _ 


del 


32v„„  dt. 


v  cG 


1  3v 


cG 


ddRdtf  dvc6  vc6  3 0R  j 


dv 


cG 


ddB 


U 


+  - 


r2  V'r 


1  3v„ 


\Vcr  }}0R 


ri  'C)Vn- 

nvt  3  0„ 


+  ■ 


2v„ 


3v 


d0D 


+  - 


ri  3vf 


\  cr 


3  0„ 


Wcr  dvr, 

Vl  d0RJ 


dv 


c0 


dt 


32v  „  dt 


v  cG 


dORdtf  dvc0 


d\ 


a'^V 


3^ 


3v 


c0 


v y 


+ 


3<T  32y 


c0 


dORdtf 


K  ^ 

dvce  ddRdtf 


+ 


32C  32^  3v 


\dvc8ddR 


-  + 


cG 


dv2c0  d  6 R  j 


dv 


cO 


dt. 


3V 


3  0R  dOR 


d2C+d£d\^  + 


dvc6  dOR 


d2C  d2C  dv 


dvceddR 


■  + 


c6 


dv;e  dOR  j 


dv 


C0 


ddv 


Differentiating  Equation  (E-8'),  namely, 


dvcff  _  1 


3/ 1 


<K 

\ddK 


+  v 


cG 


we  obtain 


i 


drldtf 


d2\ 


\  dORdtf 


-  +  - 


dv 


c0 


dt 


f  ) 


^ce  _  1 


fd\ 


dr^d  0R 


-  +  ■ 


3v 


C0 


'1  V 


del  deR  j 


Differentiating  Equation  (E-18),  namely, 


3v 

cr 

dr. 


1 


^  +  va)cotS„ 

Kdvce  3 eR 


we  obtain 


32v„.  _  1  dC  d\ 


(  ^2 


drxdtf 


ri  |  dvc6  d0Rdtf 


■  + 


\dv2ce  3 6R 


+  cot  dD 


dv 


c0 


dt , 


drlddR 


r, 


d2c  +d2Cdvc0 


V  dvc0ddR  dv]edeRj 


dvcr  ,  dC,  d2vcr  ,  3v 


■  + 


deR  dvce  de\ 


+- 


c0 


d0„ 


cot  0R  -  vc0  esc  0 R 
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d  dOR  d  ddR  d  dfd  3  dfd 

c)  6 R  dtf  c)/3h  dtf  r)  0R  dtf  d/3h  dtf 
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The  following  identities  will  be  needed*  (Appendix  G): 

cos  0R  =  sin  tp  sin  tpT  +  cos  tp  cos  tpT  cos  AD ,  (K- 1 ) 

cos  j8  sin  0R  =  cos  ^?sin  tpT  -  sin  tp cos  tpT  cos  An  ,  (K-2) 

sin  6r  sin  /3h  =  -  cos  ^sin  AD  ,  (K-3) 

cos  AD  =  -  cos  1 3  cos  Ph  -  sin  / 3  sin  /3b  cos  0R ,  (K-4) 

sin  6r  sin  /3  =  cos  (pT  sin  AD ,  (K-5) 

cos  tp  sin  p  =  -  cos  (pT  sin  /3b ,  (K-6) 

cos  tpT  cos  Ph  =  sin  ^?sin  0R  -  cos  tp cos  0R  cos  / 3 .  (K-7) 


Note  that  Equations  (K-l),  (K-5),  and  (K-7)  are  identical  to  Equations  (G-3),  (G-4),  and  (G-6), 
respectively.  It  is  clear  that  tp ,  tpT,  and  A,,  are  independent  of  each  other  so  that 

dtpT  _dtpT  _  dAD 


dtp  dAD  dtp 


■  0.  We  have  from  Equations  (G-8)  and  (G-9),  respectively, 


deR  i  ddR  .  0  dp  i  dp 

-  -cos^sinp,  ^—  =  ——r —  =  -csc0Rcos<pTcosPb. 


dAD  Q  dtf 


dAD  Q  dtf 


Differentiating  Equation  (K-l)  with  respect  to  tp  and  using  Equation  (K-2),  we  have 

•  a.  .  .  .  o  r,  d0R  0 

-sin 0R — -  =  cos^sin(Z>r  -sin^cos^  cos/tD  =  cospsin6(?  =4> — -  =  -  cos p. 
dtp  dtp 

Differentiating  Equation  (K-3)  with  respect  to  An  and  using  Equation  (K-4),  we  have 

cos  0R  sin  Ph  +  sin  0R  cos  Pb  =  -  cos  tpc os  AD 

dAD  dAD 

=4>  sin  0R  cos  Pb  =  -  cos  (p  (cos  An  +  sin  p  sin  Pb  cos  0R )  =  cos  (p  cos  p  cos  Pb 

3Ad 


dPb _ dP„  _ 


dA, 


=  CSC#RCOS^COS/? 


dtf 


=  Qcsc0r  COS  (p COS p. 


Differentiating  Equation  (K-3)  with  respect  to  (p  and  using  Equations  (K-5),  (K-6),  and  (K-7), 
we  have 

cos  0R  sin  Ph  +  sjn  q  cos  p  =  sin  sin  ; 
dtp  dtp 


Sofair,  Isaac,  K40  Training  Guide  4085.1,  A  Detailed  Derivation  of  Formulae  Arising  in  Spherical  Trigonometry, 
Naval  Surface  Weapons  Center,  Dahlgren,  VA,  May  1986  (currently  Naval  Surface  Warfare  Center,  Dahlgren 
Division). 
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sin  0R  cos  [3h  =  sin  <p  sin  An  +  cos  0R  cos  f3  sin  J3b=-^^—  (sin  (p  sin  0R  -  cos  (p cos  0R  cos  / 3 ) 

d<p  cos  (p^ 


=  sin/?  cos  (3b 


d(p 


=  esc  0Rsm(3. 


d~6> r  r/d/3  Q  Q 

— f-  =  cos  ^cos  p  =  -esc  0R  cos  (pc os  (pT  cos  p  cos  ph. 

dAD  a/ lD 

d20R  .  a  d/3  ■  ■  P  pd/3  n  ■  p  p 

- —  =  sin p—1—  =  -sin^sin p  +  cos^cos pp—  =  -esc 0R  cos (pT  sinpeos ph 

dADd(p  o/ lD  d<p 

=^>  cos  ^cos  jB  =  sin  / 3  (sin  cp  -  esc  0R  cos  (pT  cos  (3h ) 

<)(p 

=  esc  0R  sin  /?  (sin  <p  sin  0R  -  cos  cpT  cos  [3b ) 

=  cot  0R  cos  ^cos  / 3  sin  /?  (from  Equation  (K-7)) 

dj3  .  p 

=>  — —  =  cot  0R  sin  p. 

d(p 

d2 j3  d  ,  „  ,  (pa  d£?s  ■  p  d/3h  ' 

=  -cos  (pT  — —  (esc  Or  cos  ph )  =  cos  (pT  esc  Or  cos  j3b  cot  0R  -f-  +  sin  ph — A- 
oAd  dAD  (_  oAd  oAd  J 

=  esc2  Or  cos  ^cos  <pT  (cos  / 3 sin  (3h  +  sin  /? cos  /3b  cos  0R ) . 


d2j3 

dADd<p 


=  -  cos  (pT  —  (esc  Or  cos  (3h ) 
dep 

t  f  dO 

=  esc2  0R  cos  (pT  cos  Or  cos  /3h  — —  +  sin  0R  sin  y 3h 

l  d<P 

=  esc2  Or  cos  (pT  (sin  (3 sin  (3h  -  cos  (3  cos  f3b  cos  0 R 


d<p , 
)• 


It  is  desired  to  compute  ^  ,  — — — and  — — where  the  partials 

d0R  dAn  d(3h  dAn  d  Or  dAn  df3h  dAD 

with  respect  to  0 R  imply  a  constant  value  of  /3h ,  and  the  partials  with  respect  to  /3h  imply  a 
constant  value  of  0R .  Since  0R  and  [3  -  and  hence  their  derivatives  -  depend  on  (p,  <pT ,  and  An , 
these  partials  cannot  be  found  by  direct  differentiation.  In  general,  if  a  function  is  defined  by 
p(yx  (x1,x2),y2{xvx2),x1,x2)  =constant,  then 


dp_dyL+dp_dy^+dp_^0  dp  dy,  |  dp  dy2  |  dp  _() 
dy,  dxx  dy2  dv,  d.r,  d y,  dx2  dy2  dx2  dx2 
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Solving  these  two  equations  simultaneously  yields 

dp  dy2  dp  dy2  dp  dyy  dp  dy] 

dp  _  dx2  dxl  dxx  dx2  dp  _  dxt  dx2  dx2  dxx 

dyy  dyLdy^_dyLdy^'  dy2  foi  <ty2  dyx  dy2  ' 

dxx  dx2  dx2  dx]  dxx  dx2  dx2  dxx 

dO 

Putting  x,  =  (p,  x2  =  AD,  yl  =  0R,  y2  =  f3h,  and  p  =  — — ,  we  obtain 

dAD 


d20R  dph  d20R  dph 
d  dOR  _  dA2  d(p  dApdtpd/ lD  _ 
d0RdAD  d 9 R  dpb  dOR  dph 

dtp  dAn  dAD  dtp 


d  dOR 
dpb  dAD 


d20R  dOR  d20R  dOR 
dA Dd<p  dAD  dA2  d<p 
dOR  dph  dOR  dph 
dtp  dAn  dAn  dtp 


=  cos  (pT  cos  j3b. 


Similarly  by  putting  p  = 


M. 

dA/: 


,  we  obtain 


d  dp 
dOR  dAn 


d2j3  dfib  d2jB  dfib 
dA2  d(p  dADd<p  dAD 
dOR  dfih  dOR  dfib 
dtp  dAn  dAn  d<p 


esc  0R  cot  0R  cos  (pT  cos  y 3b , 


d  d/3 
dph  dAD 


d2j3  d0R  d2p  d0R 

dADd(p  dAn  dA2  d(p 

d0R  d[3h  d0R  dfib 


=  -  esc  0R  cos  (pT  sin  (5b. 


dcp  dAD  dAn  dcp 


Hence, 


d  clp 
d  9 R  dtf 


d  d0K  n  d  d0K 

- —  =  0, - - 

d  9 R  dtf  dPh  dtf 

-Q  esc  0R  cot  9 R  cos  <pT  cos  Ph , 


=  Q. cos  (pT  cos  Pb, 


d  dp 

dPi,  dtf 


-Q.  csc#Rcos^rsin/?fo. 
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